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ABSTRACT 

We present total and polarized intensity images of 15 active galactic nuclei 

obtained with the Very Long Baseline Array at 7 mm wavelength at 17 epochs 

from 1998 March to 2001 April. At some epochs the images are accompanied 

by nearly simultaneous polarization measurements at 3 mm, 1.35/0.85 mm, and 

optical wavelengths. Here we analyze the 7 mm images to define the properties 

of the jets of two radio galaxies, five BL Lac objects, and eight quasars on an- 

i/-) ' gular scales > 0.1 milliarcseconds. We determine the apparent velocities of 106 
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features in the jets. For many of the features we derive Doppler factors using 
a new method based on comparison of timescale of decline in flux density with 
the light-travel time across the emitting region. This allows us to estimate the 
Lorentz factors, T, intrinsic brightness temperatures, and viewing angles of 73 
superluminal knots, as well as the opening angle of the jet for each source. The 
Lorentz factors of the jet flows in the different blazars range from T ~ 5 to 40 
with the majority of the quasar components having T ~ 16 — 18 while the values 
in the BL Lac objects are more uniformly distributed. The brightest knots in the 
quasars have the highest apparent speeds, while the more slowly moving com- 
ponents are pronounced in the BL Lac objects. The quasars in our sample have 
similar opening angles and marginally smaller viewing angles than the BL Lacs. 
The two radio galaxies have lower Lorentz factors and wider viewing angles than 
the blazars. Opening angle and Lorentz factor are inversely proportional, as pre- 
dicted by gas-dynamic models. The brightness temperature drops more abruptly 
with distance from the core in the BL Lac objects than in the quasars and radio 
galaxies, perhaps owing to stronger magnetic fields in the former resulting in 
more severe synchrotron losses of the highest energy electrons. In nine sources 
we detect statistically meaningful deviations from ballistic motion, with the ma- 
jority of components accelerating with distance from the core. In six sources we 
identify jet features with characteristics of trailing shocks that form behind the 
primary strong perturbations in jet simulations. The apparent speeds of these 
components increase with distance from the core, suggestive of acceleration of 
the underlying jet. 

Subject headings: galaxies: active — galaxies: quasars: individual (0420-014, 
0528+134, 3C 273, 3C 279, PKS1510-089, 3C 345, CTA102, 3C 454.3) — galaxies: 
BL Lacertae objects: individual (3C 66A, OJ 287, 1803+784, 1823+568, BLLac) 
— galaxies: individual(3C 111, 3C 120) — galaxies: jet 

1. Introduction 

The jets observed emanating from the nuclei of quasars and other active galactic nuclei 
(AGN) represent the most energetic long-lived phenomenon in the universe. Although jet 
imaging is no longer the sole prerogative of radio astronomy, the exquisite resolution of Very 
Long Baseline Interferometry (VLBI) at radio wavelengths remains an unreached goal for 
sub-millimeter (sub-mm) and shorter wavelengths. It is thought that accretion onto a black 
hole drives the jet outward via magnetic forces (e.g. Meier, Koide, & Uchida 2000). Currently, 
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the most direct way to provide observational evidence for such a model is to study the linear 
polarization of jets at different frequencies (from radio to optical) along with changes in the 
innermost jet structure. The most promising objects in such investigations are blazars, flat- 
spectrum radio-loud quasars and BL Lac objects characterized by high optical polarization 
up to 46% (Mead et al. 1990; Impey & Tapia 1990), pronounced and rapid variability of 
flux, and one-sided jets with knots that move at superluminal apparent velocities, as fast 
as >30c (Jorstad et al. 2001a; Kellermann et al. 2004). According to Wills et al. (1992), 
there is a highly significant correlation between optical polarization and dominance of a 
compact radio core-jet structure, which suggests beaming of the optical flux along with the 
radio emission. Comparison of radio-loud quasars with low and high optical polarization 
shows a higher fractional polarization of the radio core for the latter (Lister & Smith 2000), 
which implies a co-spatial origin of the emission at these wavelengths. In support of this, 
the optical polarization in blazars is affected by the formation and emergence of new VLBI 
knots (Gabuzda & Sitko 1994; Gabuzda, Sitko, & Smith 1996). These results indicate 
that simultaneous multifrequency polarization monitoring, together with high resolution 
polarimetric imaging of the radio jets, provides a unique tool for identifying the location of 
the regions responsible for variability at different wavelengths and for relating the magnetic 
field geometry to the structure of the jet. 

We have obtained total and polarized intensity images of 15 AGNs with the Very Long 
Baseline Array (VLB A) at 7 mm (43 GHz) at 17 epochs over three years. The VLB A 
observations are accompanied at many epochs by nearly simultaneous (within two weeks) 
measurements of polarization at 1.35/0.85 mm (230/350 GHz) and at optical wavelengths. 
In the second half of the program simultaneous polarization observations at 3 mm were 
performed at several epochs. The main goals of the project are to relate emission regions at 
high frequencies to the parsec-scale jet structure and to investigate the strength, direction, 
and variability of the magnetic field close to the central engine. These can be achieved only 
after detailed study of the jet kinematics. This paper is devoted to an analysis of the jet 
structure and its variability associated with ejection and propagation of disturbances down 
the jet, which appear on radio maps as knots of enhanced brightness. The paper is the 
first in a series based on the entire data set collected during the project. Some results on 
individual sources have been published already by Marscher et al. (2002, 3C 120), Stirling 
et al. (2003, BL Lac), and Jorstad et al. (2004, 3C 279). Other papers will include (1) 
comparison of the polarization parameters at different frequencies with the jet structure and 
disturbances in the jet; (2) analysis of the available radio, sub-mm, optical, and X-ray light 
curves to relate the flux variability to activity in the jet; (3) results regarding stability of the 
VLBI core position based on phase-referencing observations obtained for five sources in the 
sample; and (4) structure of the core and intraday total and polarized intensity variability 
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in the parsec-scale jets. 

Although much progress has been made in studying the kinematics of jets since the 
VLBA started to operate fully in 1995 (Pearson et al. 1998; Kellermann et al. 1998; Homan 
et al. 2001; Jorstad et al. 2001a; Gomez et al. 2001; Kellermann et al. 2004), our program 
is unique regarding the number of sources observed in such detail over a rather long period 
of time. This allows us to separate the fast and slow, ballistic and curved motion in the jet 
flow, define physical parameters for both individual jet features and entire parsec-scale jets, 
and compare the results across different classes of AGNs. 

2. Sample Selection 

Our program is designed for comparison of the linear polarization at high frequencies 
(mm, sub-mm, and optical wavelengths) with the parsec scale jet structure of AGNs. This has 
defined the main criteria for selecting the sample and frequency of our VLBA observations: 

1. The sources should be bright, >0.5 Jy, and polarized, >3%, at sub-mm wavelengths. 

2. The size of the sample and brightness of sources should allow us to perform VLBA 
observations at a single epoch during 24 hours with sufficient uv-coverage to produce 
total and polarized intensity images at 43 GHz with high dynamic range. 

3. The sample should contain sources with resolved radio structure from different sub- 
classes of AGNs for which variability in the jet flow can be expected on timescales of 
months. 

4. The sources should be convenient for monitoring in the northern hemisphere and their 
coordinates should cover the whole range of right ascensions. 

Following these constraints, we have formed a sample of 15 AGNs that consists of 8 
quasars, 5 BL Lac objects, and 2 radio galaxies (3C 120 of Fanaroff- Riley type 1 and 3C 111 
of type 2) after obtaining information on the linear polarization at sub-mm wavelengths from 
Nartallo et al. (1998). The sources are listed in Table 1. 

3. Observations and Data Reduction 

We have observed the objects in our sample in four different regions of the electromag- 
netic spectrum: at 43 GHz (7 mm) with the VLBA, at 350/230 GHz (0.85/1.3 mm) with 
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the James Clerk Maxwell Telescope (JCMT, Mauna Kea, Hawaii) using SCUBA (Holland 
et al. 1999) and its polarimeter (Greaves et al. 2003), at the Steward Observatory 1.5 m 
telescope (Mt. Lemmon, Arizona) with the Two-Holer Polarimeter /Photometer (Sitko et 
al. 1985) over an effective wavelength range of ~6000-7000A, and at 86 GHz (3 mm) with 
the Berkeley-Illinois-Maryland Array (BIMA, Hat Creek, California). The observations were 
performed from 1998 March to 2001 April. The VLBA monitoring was carried out roughly 
bimonthly (17 epochs). For the JCMT and optical observations the number of epochs de- 
pends on the source, with the maximum being 11 and seven epochs and the minimum five 
and three epochs, respectively, except for the BL Lac object 1803+784, which was not ob- 
served in the optical region owing an inaccessibly high declination. The majority of the 
JCMT polarization observations were accompanied by total flux measurements. Differential 
V-band photometry (i?-band photometry in the case of 3C 273) was carried out when con- 
ditions were photometric. The BIMA polarization observations started in 2000 April and 
for many sources in the sample were performed simultaneously with the VLBA observations 
at 3-4 epochs. However, each of these epochs usually was accompanied by a series of BIMA 
observations separated by several days so that for some sources polarization measurements 
at 3 mm were obtained at 20 epochs. The polarization data at 3 mm include all four Stokes 
parameters. In this paper the optical, JCMT, and BIMA polarization results are shown only 
graphically to illustrate the data collected. Detailed descriptions of the data reduction and 
tables of measurements will be presented in a later paper. 



3.1. Radio Observations 



The 43 GHz observations were carried out with the VLBA recording system using eight 
8 MHz wide channels, each in right and left circular polarization, with 15-20 scans of 3-5 
minute duration for each object. All 10 antennas were used for each source except at epochs 
affected by weather or receiver failure. Table 2 lists the antennas operated at each epoch and 
typical parameters of the synthesized beam (those of the OJ 287 observations) for uniform 
weighting (used for imaging all sources) to illustrate the consistency of the liv-coverage over 
epochs. Initial correlation was carried out at the National Radio Astronomy Observatory 
(NRAO) Array Operations Center in Socorro, NM. Subsequent calibration was performed 
with the Astronomical Image Processing System (AIPS) software supplied by NRAO, while 
images were made with the Caltech software Difmap (Shepherd 1997). The calibration 
included application of the nominal antenna-based gain curves and system temperatures, 
and correction for sky opacity, followed by iterative imaging plus phase and amplitude self- 
calibration. For each epoch we calculated the total flux density in the images of sources that 
are known to have very weak emission outside the angular size range of the VLBA images 
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(0420-014, 0528+134, OJ 287, and BL Lac). These values were compared with total flux 
densities obtained by interpolating in time the measurements of the monitoring program at 
37 GHz at the Metsahovi Research Station, Finland (Terasranta et al. 2004). The comparison 
produced the flux density correction factors, / amp , given in Table 2. The factors were applied 
for the final adjustment of the flux density scale in the images. We have constructed light 
curves of the VLBI core for each source to confirm the absence of correlation between their 
behavior, indicating no evidence for residual amplitude calibration errors. The light curves of 
the 15 sources peak at 9 different epochs and the correlation coefficients are spread between 
—0.4 and +0.3, consistent with no systematic calibration errors in the flux density scaling. 

We performed a cross-hand fringe fit using a scan of 3C 279 averaged over all baselines 
(the data were corrected for parallactic angle rotation of phases in advance). The resulting 
right-left phase and rate delay corrections were processed with the AIPS task POLSN and 
applied to the full data set. After preliminary images were produced, the images were used 
in AIPS task CALIB to self-calibrate the phases of all sources; this includes the removal of 
residual R-L phase differences owing to calibration errors. The instrumental polarization "D- 
terms" were determined via the method of Roberts, Wardle, & Brown (1994) and Leppanen, 
Zensus, & Diamond (1995). A final set of D-terms at each epoch was obtained by averaging 
of the solutions for those sources in the sample for which there is the best agreement between 
the D-terms (usually, 3C 111, OJ 287, 3C 345, and BL Lac). The electric vector position angle 
(EVPA) calibration was obtained by different methods: comparison between the Very Large 
Array (VLA) and VLBA integrated EVPAs at quasi- simultaneous epochs, D-terms method 
(see below), and using EVPA-stable features in the images of the jets of 3C 279, OJ 287, and 
CTA102. Over the period 1998 March - 1999 February we used results of the VLA observa- 
tions performed for 0420-014, 3C 120, OJ 287, 3C 279, BL Lac, and 3C 454.3 by Gomez et 
al. (2002). Later we obtained polarization measurements for 0420—014, 0528+134, OJ 287, 
and 1803+784 at two epochs, 2000 April and 2000 July, with the VLA C and CnD config- 
uration, respectively. At epochs where VLA data were not available we used the D-terms 
method, which is based on the assumption that the instrumental polarization parameters 
change slowly with time (Leppanen, Zensus, & Diamond 1995). A detailed description of the 
method is given in Gomez et al. (2002). In addition, later we checked our calibration for these 
epochs using the NRAO data base at http:/ /www. vla.nrao.edu/astro/calib/polar/ for 
sources 0420-014, 0528+134, OJ287, and BL Lac when the NRAO data were available at 
epochs close to the VLBA dates. The results are completely consistent with our calibration. 
The EVPA-stable feature in 3C 279 is the well-known superluminal component C4 (e.g., 
Homan et al. 2003), which underwent a change in EVPA in the beginning of 1998 (Jorstad 
et al. 2004) but maintained this EVPA until the end of our program. The EVPA-stable 
features in OJ 287 and CTA 102 are quasi-stationary components at ~1 mas and ~2 mas 
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from the core, respectively. However, these both became very weak in the second half of 2000 
in the 43 GHz images. The final EVPA calibration is a result of the best agreement among 
the different methods. The accuracy of the calibrated EVPA measurements, o"(EVPA), is 
indicated in Table 2. 

We have combined the images obtained for every source in a sequence of total intensity 
maps convolved with the same beam, corresponding to the average beam during the epochs 
when all 10 antennas were in operation. The contours are in terms of the global peak 
of the maximum total intensity observed over all epochs. The sequences are presented in 
Figures 1-15, which also show the JCMT, BIMA, and optical polarization measurements. All 
images are oriented with north toward the top and east toward the left. The scale in mas is 
indicated along one of the axes, while the epoch of each image is given along the other axis. 
The peaks of the total, I pca k, and polarized, /p eak , intensity, parameters of the average beam, 
the total intensity of the lowest contour, J min , and the lowest level of the polarized intensity, 
7^ in , of the combined maps are indicated in Table 1. The polarization at 7 mm is shown 
inside the total intensity images by line segments, which are plotted if the local polarized 
intensity exceeds i^m- The segments are oriented in the local direction of the polarization 
and their length is proportional to the local polarized intensity, with the maximum length 
corresponding to the global peak of the polarized intensity over all epochs, /p eak . Table 2 
marks epochs when the BIMA, JCMT, and optical observations were performed within two 
weeks of the corresponding VLBA epoch. 



3.2. Model Fitting of the VLBA Images 

We employed the task MODELFIT in Difmap to represent each total intensity image 
(Stokes parameter /) as a sequence of circular Gaussian components that are characterized 
by flux density, size, and position relative to the map center. Initially, point-like compo- 
nents are used to obtain an image with 100:1 dynamic range. When a set of point-like 
components was found, a final group of 100 iterations was executed with all parameters 
of all components allowed to vary to define their properties. In addition, because we have 
roughly bimonthly observations, we repeated the model fitting with the components from 
the previous epoch used as the initial model (except for the first epoch, 1998 March). This 
improves identification of components over epochs and provides an estimate of the accuracy 
of the parameters by comparison of the outcomes obtained with different initial models. For 
several sources (3C 66A, OJ 287, 3C 279, PKS 1510-089, and 1803+784) at four epochs 
(1998 March 1998, 1999 October, 2000 July 2000, and 2001 January) we derived estimates 
of the la uncertainties of the best-fit models using the method described by Tzioumis et al. 
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(1989) and realized in the package Difwrap (Lovell 2000). The uncertainties depend signifi- 
cantly on the brightness and size of components: 1) for bright (flux of knot > 100 rms noise 
level) and compact (size <0.1 mas) features, uncertainties in flux density ~1%, in position 
~0.01 mas, and in size < 1%; 2) for the majority of components with size 0.1-0.3 mas and 
flux >50 mJy, uncertainties in the flux density ~3%, in position ~l/5 of the beam size, and 
in size ~5%; 3) for diffuse components, uncertainties in flux density ~10%, in position are 
comparable with half the size of the knot, and in size ~10%; 4) for components weaker than 
50 mJy, uncertainties in flux density ~50% and positional uncertainties correspond to the 
size of the average beam (see Table 1); however, for the majority of weak components these 
uncertainties were later modified as discussed in §4 (Errors of Polynomial Parameters). 

We determined the parameters of polarized components by applying the task MOD- 
ELFIT to the uv-data in Stokes parameters Q and U separately. For this purpose only 
point-like components were used in the models. Because the Q and U components can be 
either positive or negative, the procedure needs to be carried out carefully by checking with 
the polarized intensity image to avoid the inclusion of false features generated by the imaging 
procedure. The derived positional parameters of the Q and U components were compared 
with the parameters of / components. We consider a Q and/or U component to be associ- 
ated with an I component if they are co-spatial to within the uncertainties. If both polarized 
components are present, the position of the polarized component is defined by the average 
location, weighted by the absolute values of the Q and U fluxes. There are few cases when a 
significant polarized component does not have an / counterpart. Uncertainties in the values 
of the parameters of polarized components are difficult to define. Estimates of the accuracy 
of polarization parameters in jet features at 15 and 22 GHz are discussed by Homan et al. 
(2002a). For components that are bright in the total intensity and highly polarized (percent 
polarization >5%), we find that the uncertainties in fractional polarization ~ 1% and in 
polarization angle ~ 5°. 

A table of parameters of jet features for each source over all epochs can be found at the 
web-site www.bu.edu/blazars/multi.html. The columns of these tables are as follows: 
1 - epoch, 2 - flux density in Jy, 3 - relative right ascension in mas, 4 - relative declination 
in mas, 5 - distance from the core in mas, 6 - position angle relative to the core in degrees, 7 
- angular size in mas, 8 - polarized flux density of a polarized component associated with a 
total intensity component, in Jy, 9 - distance of polarized component from the core in mas, 
10 - position angle of polarized component relative to the core in degrees, 11 - EVPA in 
degrees. 
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4. Technique to Measure Motion in the Jets 



The data allow us to follow the evolution of major features of the jets that are observed 
in the total and polarized intensity VLB A images of the 15 AGNs at 17 epochs over 3 years. 
For each image we identify a component, AO, as the VLBI core. For all sources AO is located 
at one end of the jet, although it is not always the brightest feature. The position of the core 
in right ascension x and declination y is defined as x = and y — 0. In our analysis we assume 
that the core is stationary over all epochs. Locations of other features are determined relative 
to the core. Components are categorized as follows: knots A (other than AO) are features 
that either are stationary (the proper motion is less than or equal to its uncertainty), undergo 
reverse motion (toward the core), or move at subluminal apparent speeds. Knots B, C, and D 
are superluminal features, where B components are the fastest knots in the jet, D components 
are the farthest knots from the core, ejected before our monitoring period started, and the 
remainder are labeled as C. The designation of components is different for 3C 120, 3C 279, 
3C 345, and BL Lac, where the naming corresponds to Gomez et al. (2001), Wehrle et 
al. (2001), Ros et al. (2000), and Stirling et al. (2003), respectively. Each component is 
characterized by the following parameters: total S and polarized S p flux density; positions x 
(RA), y (Dec), and R, where R = \J x 2 + y 2 ; position angle = t&n~ 1 (x/y); size (FWHM) 
of component a; and electric vector position angle, EVPA. To define the temporal evolution 
of the jet and determine the apparent velocities of the jet flow we perform the following 
steps. 

1. Identification of components at different epochs. The identification of a component is 
based on comparison of the parameters during the epochs when it is visible on the images. 
We assume that the same component has similar total, S, and polarized, S p , fluxes, position 
angles 9 and EVPA, and size, a, at successive epochs. However, some components evolve 
dramatically even over a two-month interval between epochs, splitting into one or more 
subcomponents or merging with other features of the jet. In this paper we analyze those 
jet features whose identification is supported by the similarity of a number of parameters at 
different epochs. 

2. Fitting Various Polynomials. The thorough sequences of images allow us to search 
for acceleration/deceleration of the jet flow and non-ballistic projected trajectories. We fit 
the x,y positions of a component over N epochs by different polynomials of order /: 



where U is the epoch of observation, i=l,...,N, and t mi d=(ti + t N )/2. We use the program 
LSQPOL of the FORTRAN version of the package DATAN (Brandt 1999) to find the optimal 




(2) 
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polynomials of order /, where / runs from to 4. The upper limit, I — 4, is a consequence 
of the maximum number of epochs of observation of a knot, equal to 17. The program 
provides the value Mi, which is the goodness-of-fit by a polynomial of order I. For each 
order, we perform a x 2 test to determine the polynomial that best fits the data. We choose 
a polynomial of the lowest order for which Mi < M x 2, where M x 2 is the value of the x 2 
distribution corresponding to significance level C=0.05 for / = N — I — 1 degrees of freedom 
(Bowker & Lieberman 1972). 

Examples of selection of polynomials that fit the data are given in Tables 3 and 4 for 
component CI identified in the jet of 3C 66A at 16 epochs and component Bl, also seen 
at 16 epochs, in 3C 273. These tables show parameters of the best-fit polynomials of order 
from to 4, the corresponding values of the goodness of fit to the data M, and the goodness 
required by the x 2 -test. Table 3 shows that 1=1 satisfies the x 2 test for coordinates x and 
y in the case of CI for 3C 66A. For Bl, l=A is required to match motion in RA, while a 
second-order polynomial adequately describes the data for declination. 

In a few cases the x 2 test is not satisfied even by polynomials of order 4. An increase in 
the order of the polynomial does not improve the situation since, given the limited number 
of observations, the number of degrees of freedom becomes too small. Such cases are ap- 
proximated by a straight line with special consideration for uncertainties in the x,y values 
(see below). These cases are marked in Table 5 by "1*" under the polynomial order. 

3. Uncertainties in Polynomial Parameters. After the best-fit polynomial is found, 
we run the subroutine LSQASN of the package DATAN to estimate errors in the derived 
parameters. This method is valid for a normal distribution of the unknowns and when 
the true value lies with probability W within the confidence region around the estimated 
value (W=0.95 is used). However, the errors depend significantly on the uncertainties in 
the individual observations that are obtained from the model fitting (see §3). This allows us 
to lessen the weight of data at epochs with bad weather, a failure of one or more antennas, 
or when a given component was very weak or diffuse. For component motion that cannot 
be represented well by polynomials I <4 (see above), we approximate the motion by a first- 
order polynomial and determine uncertainties in the x,y position using the program LSQPOL 
and the method suggested by Homan et al. (2001). We set the uncertainty for each data 
point equal to the average beam size and calculate the optimal polynomial of first order 
corresponding to a preliminary x 2 value. Taking this preliminary x 2 value, we uniformly 
re-scale the uncertainties in the data points such that x 2 would correspond to the x 2 value 
for (=0.05 and / = N - 2. 

4. Calculation of Proper Motion, Acceleration, and Ejection Time. We define the 
average proper motion as a vector (< /i >,< $ >), where < /i > represents the mean 



- 11 - 



angular speed of motion and < $ > gives the average direction of motion. The average 
values are derived as follows: < /i >= ^ < fi x > 2 + < fi y > 2 and < $ >= tan -1 ^^ ), 

where < \i x >= f** xdt/(tN — h) and < fi y >= f** ydt/(tN — ti). In the case of first or 
second-order polynomials, < fi x >—a\ and < \i y >=b±. For a polynomial of order higher than 
unity we compute an average vector of acceleration, (fi»,fi±), where fi\\ is along the direction 
of the average velocity, < $ >. In the case of a second-order polynomial fi x = 2 a 2 and 
fly = 2 6 2 - For a higher order polynomial, fi x = f* N xdt/(tN — ti) and fi y = f^ N ydt/ (t^ — ti). 



5. Jet Velocities 

We find superluminal apparent speeds for 19 of 22 components in the radio galaxies, 19 
of 31 knots in the BL Lac objects, and 46 of 53 knots in the quasars. Table 5 lists the average 
apparent speed, < f3 app >, calculated in units of the speed of light, c, using the average proper 
motion, < fj, >. An inhomogeneous Friedmann-Lemaitre- Robertson- Walker cosmology, with 
Vl m = 0.3, Ha = 0.7, and Hubble constant H o =70 km s _1 Mpc -1 (Kantowski, Kao, & 
Thomas 2000), is adopted for the calculations. We derive the uncertainties in the average 
proper motion, acceleration, and apparent speed from the uncertainties in the polynomial 
coefficients. Table 5 contains also the number of epochs, N, at which the component has been 
observed; the average total flux, < S >, in Jy; the average distance from the core, < R >, 
in mas; the average position angle of the component, < >, in degrees, with uncertainty 
corresponding to the scatter across the epochs; the average direction of the velocity vector, 
< $ >, in degrees; and the ejection time (epoch of zero separation), T Q . The ejection time 
is the extrapolated time of coincidence of the position of a moving knot with the core in the 
VLBA images. T Q is the average of t x0 and t yo weighted by their uncertainties, where t xo and 
t yo are roots of the best-fit polynomials. In the case of a polynomial of order 3 and higher, 
the roots are computed by successive iterations with an accuracy of 10~ 5 . 

Table 6 gives the parameters of acceleration for components in Table 5 having best-fit 
polynomials of order >2 (some columns repeat information from Table 5 for convenience). 
Figure 16 shows the positions of the jet components relative to the core for both the to- 
tal (open circles) and polarized (filled circles) intensity images of each source. The solid 
lines/curves indicate the best polynomial fit to the data. Components presented in Figure 
16 are marked at the first epoch in Figures 1-15 if they are detected then or at the later epoch 
at which they appear in the jet. Figure 17 displays an image from one particular epoch for 
each source when the most prominent jet features are seen. Table 7 lists the parameters of 
the maps shown in Figure 17 (rms(I) and rms(Ip) are the root mean squares of the residual 
total and polarized intensity on the images, respectively). The trajectories of all components 
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with non-ballistic motion are plotted in Figure 18. Appendix A contains a description of our 
results for each individual object in the sample. 

6. Apparent Speed as a Probe of Parsec-Scale Jets 

Our intensive, prolonged high-frequency VLBI monitoring reveals complex, changing 
structure of parsec scale jets, including significant variations in the apparent speed of knots 
in the jet, /3 apP , within individual objects. The apparent speed is determined by the intrinsic 
velocity, f3, and the angle between the trajectory and the line of sight, 6 : 

l3 ap p = 13 sin O (1 - (3 cos O ) _1 . (3) 

The intrinsic velocity defines the Lorentz factor of the knot, T = 1/a/1 — (3 2 , where (3 is 
in units of the speed of light. Large scatter in the apparent speed within a source can be 
caused by (1) different patterns of jet components, such as "blobs" of energetic plasma or 
forward, reverse, or stationary shocks; (2) variable Lorentz factor of the jet flow; or (3) 
different trajectories of components in the jet. The latter can result from a change in the 
direction of the jet or in the paths of components if each given knot does not fill the entire 
cross-section of the jet. 

The observed flux density of a superluminal component is boosted in the direction of 
the observer by a factor S 3+a , where 

5= [r(l-/3 cos Go)]" 1 (4) 

is the Doppler factor and a is the spectral index (S u oc u~ a ) of the knot. The flux of a 
component decreases with distance from the core as the result of radiative energy losses and 
expansion ("adiabatic" cooling). The emergence of a new component is associated with one 
or more flares in the radio light curves (e.g., Savolainen et al. 2002). The flare usually has a 
sharp peak and nearly symmetric exponential rise and decay at high frequencies (Terasranta 
& Valtaoja 1994). In the case of the shock-in-jet model for the appearance of superluminal 
knots, variations in flux, spectral energy distribution, and polarization (Marscher & Gear 
1985; Hughes, Aller, & Aller 1985), the symmetric light curves of flares suggest that the 
flux variability is controlled by light travel delays across the shocked region (see Sokolov, 
Marscher, & McHardy 2004). This assumption allows us to calculate the Doppler factor for 
each superluminal component using its flux density variability and size measured from the 
VLBA data. We can then use the derived Doppler factors to study the physical parameters 
of the jets in our sample. 
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6.1. Physical Parameters of the Jet Components 

Figure 19 shows light curves of all superluminal components listed in Table 5 (except 
knots classified as trailing features, see §6.6). Each light curve is normalized by the cor- 
responding average flux indicated in Table 5. We define the timescale of the variability 
for each superluminal component as At var = dt /ln(S max / S min ) (Burbidge, Jones, & O'Dell 
1974), where S max and S m i n are the measured maximum and minimum flux densities, respec- 
tively, and dt is the time in years between and S ni i n . The variability Doppler factors 
are derived as 

^var = T~ 7Z j T, (5) 
C At var (1 + Z) 

where D is the luminosity distance, s is the angular size of the component, equal to 1.6a for 
a Gaussian with FWHM=a measured at the epoch of maximum flux if the true geometry 
is similar to a uniform face-on disk. In the case of a point-like knot, we adopt a=0.1 mas, 
which yields an upper limit to 5 var . This definition of the Doppler factor assumes that the 
variability timescale corresponds to the light-travel time across the knot. This will be true 
if the radiative cooling time is shorter than the light crossing time, which in turn is shorter 
than the timescale for cooling by adiabatic expansion. We can roughly verify this assumption 
by using the relation between flux density and size of a shocked region derived in Marscher & 
Gear (1985) during the "adiabatic" stage. We calculate the timescale, At a , of the variability 
in size for each non-point-like component as At a = dt / ln(a max / a niin ) (we use a max = as min an d 
a min = a s max ^ a s min > a s max ; which is valid for the majority of components). Figure 20 plots 
the size variability timescale versus the flux variability timescale. The straight line indicates 
the expected relation between At a and At var for adiabatic losses for optically thin shocked 
gas with a = 0.7 (Marscher & Gear 1985). Figure 20 shows that the majority of components 
have shorter flux variability timescales than those predicted for adiabatic expansion. This 
implies that at high radio frequencies the decay in flux is driven by radiative losses and, 
therefore, equation 5 should apply to our data. 

A widely used method to estimate the Doppler factor from VLBI data assumes that the 
highest apparent speed detected in a source defines the lowest possible Lorentz factor of the 
jet, T > 1^/1 + /3 app . The viewing angle is then taken to be o < sin _1 (l//3 a pp); this leads 
t° ^/3a PP ~ Aipp- in Figure 21 Doppler factors thus derived are plotted versus the Doppler 
factor computed from the flux variability and size of components having the highest apparent 
speed in the jets. Figure 21 demonstrates that there is reasonable agreement between the 
Doppler factors estimated by these two different methods. The best least-square linear fit 
to the dependence is 5p &pp = (0.72 ± 0.15)o~ var , consistent with S var being the true value and 
5p being a lower limit taking in account that 5p is estimated for the lowest possible 
Lorentz factor. We conclude that the combination of variability of flux and measurement of 
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angular sizes provides a new, robust method for deriving Doppler factors from well-sampled 
sequences of VLBI data. 

Tables 8, 9, and 10 give estimates of the Lorentz factor, viewing angle, Doppler factor, 
observed brightness temperature, T bobs , and intrinsic brightness temperature, T bint , for 43 
knots in the quasars, 19 knots in the BL Lacs, and 15 knots in the radio galaxies. The Lorentz 
factor and viewing angle are solutions of the system that combines equations (3) and (4) 
for a knot having measured apparent speed with the Doppler factor derived via equation 
(5). The observed brightness temperature is computed based on the VLBI measurements as 
2b,obs = 7.5 x 10 s S m£LX /s 2 K, where .Smax is the maximum flux of the component in Jy and s 
is its measured angular size in mas at the epoch of maximum flux (see above). We apply the 
Doppler factors derived with equation (5) to estimate the intrinsic brightness temperature, 
T bli „t = T b , obs (1 + zf- 7 /5 l -\ where we adopt a = 0.7. 

6.2. Opening Angle of the Jet 

We assume that jets have conical structure and constant angle, 9, between the jet axis 
and surface of the cone that contains the entire region of emission. Therefore, 9 is the actual 
half opening angle of the jet. This assumption is likely correct for the section of the jet within 
1-2 mas from the core, while at larger distances (>10 mas) the jet could trace out a helical 
structure covering a wide range of position angles even if the body of the jet remains narrow 
(e.g., Lister et al. 2003). We estimate the projected half opening angle, 9 P , for each source 
using the ratio between apparent transverse size, s t , of the jet and apparent longitudinal 
distance, s\, of components: 9 p = tan -1 ip, where is the slope of the best linear fit to the 
relation between s t and s\ as defined at the position of each component that is brighter than 
1% of the peak intensity. The values of s t and s\ are calculated as follows: s\ = R, where R 
is the observed separation of the component from the core; and s t = R sin (|6j Ct — 0|) +a/2, 
where 0j et is the projected direction of the jet as defined by the mean position angle of all 
sufficiently bright components over all epochs, is the position angle of the component, 
and a is the size of the component. We first plot the values of s t against s\ to estimate a 
preliminary linear dependence. The preliminary relationship is used to remove points that 
deviate from the dependence by more than 3<r of the linear model. Then the value of ip is 
obtained by minimizing \ 2 of the adjusted data. In Figure 22 all pairs of (s t , s{) are plotted 
and the best linear fit of the adjusted data is presented for each source. Figure 22 shows 
that a jet model with constant opening angle provides a good approximation to the inner 
jet, although the plots of two BL Lac objects (3C 66A and 1803+784) display an increase in 
9 P beyond 2 and 1 mas, respectively, which might be a general feature of many BL Lac jets. 
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In the radio galaxy 3C 111 and quasar 3C 345 a decrease in 9 P is observed, and is most likely 
caused by the weakness and diffuse nature of components at large distance from the core, 
such that our images only contain portions of these features. Table 11 gives the parameters 
©jet; 9 P , < @o >, 9, < T >, and < 5 > for each source. Parameters < T >, < 5 >, and 
< ©o > are weighted averages of the values for individual components listed in Tables 8, 
9, and 10, with weights inversely proportional to the uncertainty in apparent speed. The 
intrinsic half opening angle, 9, is estimated as 9 = 9 P sin < O >, where < O > is the angle 
between the jet axis and the line of sight. 



6.3. Jet Parameters 

We have constructed the dependence between the Doppler beaming factor, <5 var , and 
apparent speed (Fig. 23) to compare derived parameters with results from the 2-cm VLBA 
survey (Kellermann et al. 2004). These authors determined Doppler factors using the vari- 
ability method of Lahteenmaki & Valtaoja (1999), with an intrinsic brightness temperature 
of 2 x 10 10 K that they assume corresponds to the region of the VLBI core where the main 
variability at radio frequencies occurs (Savolainen et al. 2002). In Figure 23 the majority of 
points follow the expectation that f3 app < 5 VSLT and lie inside the u l/fl app cone". Moreover, 
despite our sample consisting of blazars or blazar-like sources for which a higher Lorentz 
factor is expected, the upper limit to f3 app corresponding to V = 25 (Kellermann et al. 2004) 
applies to our sample as well (dotted curve in Fig. 23). Figure 24 shows the distribution 
of Lorentz factors of superluminal components. There is a significant scatter in T for the 
quasars and BL Lacs, and the Lorentz factors of the radio galaxies are lower than for the 
other classes. The distributions, in general, agree with the distribution of Lorentz factors 
found in the 2-cm survey. 

Monte-Carlo simulations for flux limited samples of radio sources predict a connection 
- although not direct correlation - between the Lorentz factor and viewing angle of the jets: 
the higher the Lorentz factor, the smaller is the viewing angle (Lister & Marscher 1997). 
Our sample of the blazars and blazar-like sources shows a significant correlation between 1/T 
and O (coefficient of correlation 0.83, Fig. 25); this supports the above prediction. Small 
intrinsic bends in jets oriented very close to the line of sight should be greatly amplified 
in projection on the sky (Readhead et al. 1978). The expected morphology has not been 
confirmed by VLBI surveys. For example, 7-ray blazars, which exhibit the highest apparent 
speeds (Jorstad et al. 2001a; Kellermann et al. 2004), do not possess more pronounced 
bends than sources not yet detected in 7-rays (Kellermann et al. 1998; Jorstad et al. 2001a). 
Estimates of the viewing angle obtained for our sample allow us to test the relation between 
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morphology and viewing angle of the jets more directly. We have determined the average 
projected position angle of the jet, 0j et , and its standard deviation, cr(0j et ), for each source 
(see Table 11). In sources with jet direction close to the line of sight, a small intrinsic 
change of the component trajectory should result in a significant bend in the projected 
path, leading to a large scatter in 0j et and yielding a high value of a(Oj e t)- A similar test 
is to compare the projected opening angle with the viewing angle of the jets: jets viewed 
very close to the line of sight should have a wide projected opening angle on VLBI maps. 
Figure 26 presents both plots. The left panel shows the relationship between 6 and cr(9j e t)- 
Although there is a general decrease in <x(Oj Ct ) with increasing angle between the jet axis 
and line of sight (coefficient of correlation —0.58), the relationship is very weak or absent for 
viewing angles < 5° (coefficient of correlation —0.18). The right panel shows that there is 
no connection between the projected opening angle and viewing angle of the jet (coefficient 
of correlation —0.16). This result implies that amplification of the projected size by a small 
angle between the jet axis and line of sight is partly (or sometimes completely) canceled by 
the proportionality between opening angle and 1/T (see below). Possible reasons behind the 
lack of correlation between degree of apparent bending and small viewing angle include: (1) 
an inverse relationship between bulk Lorentz factor and degree of intrinsic bending, given 
that the momentum of the jet is proportional to T, (2) broader opening angles in slower jets, 
with filamentary structure appearing similar to bending of the axis, and (3) the resolution 
available transverse to the jet not being sufficient to detect bending in many high-r objects. 

The Lorentz factor plays a significant role in determining the jet geometry. According 
to standard models of relativistic jets, the opening angle of the jet should be inversely pro- 
portional to the Lorentz factor (e.g., Blandford & Konigl 1979). In Figure 27, the estimated 
half opening angles are plotted versus the derived Lorentz factors. There is an obvious de- 
crease of 9 toward higher values of V. According to a \ 2 test the observed dependence can be 
described by the relation 6 ~ p/Y rad, where p = 0.17±0.08 (solid line in Fig. 27) at a 2.5% 
level of significance (x 2 =6.875, f=2). The current models for the formation of relativistic 
jets that employ confinement of a jet by magnetic forces (e.g. Meier, Koide, & Uchida 2000; 
Vlahakis & Konigl 2004) do not yet specify the opening angle of the jet as a function of the 
model parameters. Such an expression, however, has been derived based on gas dynamics for 
a relativistic jet confined by pressure equilibrium with its surroundings (Daly & Marscher 
1988). In this model the opening angle depends on the Lorentz factor of the flow and the 
ratio of the external pressure, P ex t, to the initial pressure, P a , of the plasma in the core 
region, £ = ^P ext /P Q . Using equations (18), (19a), and (19b) in Daly & Marscher (1988), 
we calculate dependences of the opening angle on Lorentz factor for different values of £ 
(dotted curves in Fig. 27). The best fit of the observed dependence (5% level of significance) 
coincides with the model in which £ = 0.6 (x 2 =4.586, f=2), corresponding to a pressure ratio 
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in the core region P ext /P 1/3. 

Figure 28 presents the distribution of intrinsic brightness temperature of the jet com- 
ponents. Note that T b int is measured at distances ranging from several parsecs to a few 
kiloparsecs from the core (see Fig. 28, right panel), where the knots are generally optically 
thin. The distribution peaks at T h - mt ~ 2 x 10 9 K for the quasar knots and T bi mt ~ 6 x 10 7 K 
for the BL Lac knots, while for the radio galaxies the temperatures are evenly distributed 
between these values. The maxima should be adjusted to slightly higher temperatures since 
for 30% of the quasar and BL Lac components we have obtained only upper limits to the 
Doppler factors. Comparison of these brightness temperatures with the intrinsic equipar- 
tion brightness temperature of the optically thick part of the jet (7b,int ~ 2 — 5 x 10 10 K), 
thought to depend weakly on source parameters (Readhead 1994; Lahteenmaki & Valtaoja 
1999), implies a faster drop of the intrinsic temperature from the compact core region to the 
more extended structure in the BL Lac jets. One possible explanation for a lower intrinsic 
brightness temperature is the presence of a stronger magnetic field (Readhead 1994). This 
would lead to more severe radiative losses and weaker jet components relative to the core 
in the BL Lac objects, and perhaps a lower intrinsic brightness temperature in the cores as 
well. 

In Figure 29 we plot the average Doppler factor vs. average viewing angle (left panel), 
and the average viewing angle vs. the intrinsic half opening angle of the jet for each source 
(right panel). The bold crosses indicate the average position of the quasars, BL Lac ob- 
jects, and radio galaxies, with the size of each cross equal to the la uncertainty of the 
parameters. Although each subclass has a large scatter around the average values and sta- 
tistically the difference in the parameters of the quasars and BL Lac objects is negligible, the 
points form a continuous sequence on a 3-D plot of 5, , 9 with: quasars (< 5 >=23±11, 
< O >=2.6°±1.9°, #=0.5°±0.3°) — ► BL Lacs (< 5 >=13.5±6.7, < O >=4.4°±3.0°, 
fl=0.6°±0.4°) — ► radio galaxies (< 5 >=2.8±0.9, < O >=19.5°±3.2°, #=3.2°±0.5°). 

6.4. Accelerating/Decelerating Flow 

According to the x 2 test, one component in the radio galaxy 3C 120, five components 
in two BL Lac objects, and 13 components in six quasars exhibit a statistically signifi- 
cant change in the proper motion with time corresponding to acceleration/deceleration with 
distance from the core. Therefore, in 9 out of 15 sources the apparent speed of individual 
components varies. Moreover, although 86 knots are classified as moving ballistically, ~ 38% 
have trajectories that we are unable to fit well by a polynomial of any order, causing us to 
suspect non-ballistic motion for these as well. Figure 30 shows the evolution of apparent 
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velocity, including direction along the jet, derived from the best polynomial approximation 
of component positions across epochs for knots with a detected change in apparent speed. 
We have analyzed the results to search for a general trend in the variation of apparent speed 
with distance from the core. For each component, we compare the instantaneous apparent 
speed at different distances from the core with the average (indicated in Table 5), and assign 
a value representing the apparent speed at each angular distance of +1 if the velocity is 
higher than the average, —1 if it is lower than the average, and if it equals the average. 
Then we construct the distribution of such changes along the jet using the derived viewing 
angle for each component (Tables 8-10), which allows deprojection of observed distances 
from the core. The result is shown in Figure 31 where the distribution above the x-axis 
represents accelerating components and the distribution under the x-axis indicates deceler- 
ating components at corresponding deprojected distances from the core. Figure 31 reveals 
that at distances larger than ~5 pc from the core an increase of apparent speed is more 
common. However, there are some jet components that undergo alternating periods of both 
acceleration and deceleration. The latter could be a signature of helical motion (e.g., Denn, 
Mutel, & Marscher 2000; Tateyama & Kingham 2004) or pinch instabilities that cause the 
jet cross-section, pressure, and Lorentz factor to oscillate with distance from the core (e.g., 
Gomez et al. 1997). 

The cause of the accelerations could be bending combined with selection of objects with 
high Doppler factors such that the mean angle to the line of sight of the region near the 
core is less than optimal for superluminal motion. Statistically, such jets are more likely to 
bend away from the line of sight, thus increasing their apparent speeds with distance from 
the core. Alternatively, the acceleration could be physical, caused by considerably higher 
energy density in relativistic particles than in rest mass (Daly & Marscher 1988) or magnetic 
acceleration (Vlahakis & Konigl 2004). 

6.5. Forward and Reverse Shocks 

Our data reveal the coexistence of both very fast and much slower (but moving) knots in 
the jets of the majority of the sources in the sample (see Fig. 16). Diversity in the apparent 
speeds of jet features might reflect intrinsic variations in the pattern speed of disturbances in 
the jet flow. In the shock wave model for interpretation of the radio light curves and features 
propagating down a relativistic jet, a disturbance such as an increase in velocity or energy 
flux in a jet can create both forward and reverse shocks (e.g, Hughes, Aller, & Aller 1991). 
Both move away from the central engine, but the reverse shock has a lower velocity than 
does the forward shock. In this context fast jet features can be associated with the forward 
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shock and slow moving knots with the reverse shock. Analysis of the brightness of these 
features might reveal a prevalence of different types of shock waves in the jets for different 
classes of AGNs. 

For each knot in Table 5 (except knots classified as trailing features; see §6.6) we have 
computed parameter F re i = S max / < Sao >, which characterizes the flux of a knot relative 
to the core, where is the maximum observed flux density of the knot and < Sao > is 
the average flux density of the core over epochs. Figure 32 shows the distributions of the 
derived values of this parameter for fast (f3 app > 3c) and slow (/9 a pp < 3c) moving knots in the 
quasars, BL Lac objects, and radio galaxies. The separation into fast and slow features is not 
strictly defined, and the distributions do not change significantly if the dividing velocity is 
in the range from 2c to 4c. The distributions of F rel of fast knots in the quasars and BL Lacs 
are different at 99.5% level of confidence (/=5) according to the x 2 test. The brightness of 
~50% of the fast knots in the quasars is comparable to the brightness of the core (F re i >0.5), 
while the brightness of fast knots in the BL Lac objects never exceeds half the brightness 
of the core, and the distribution peaks at F rc i <l/4. The distributions of F re i in the radio 
galaxies cannot be classified due to the small number of objects. 

The distributions for slow moving knots indicate the presence of two populations: bright 
(F re i >1) and faint (F Te \ <0.5 for the quasars and radio galaxies and F m \ <1 for the BL Lacs). 
The population of bright slow knots consists of jet features with subluminal apparent speeds 
that most likely represent stationary shocks in the jets. This population is most prominent 
in the quasars. The population of faint slow knots is prominent in the BL Lac objects, 
including 39% of all identified jet features, while in the quasars such features comprise only 
7% of detected knots. About 70% of faint slow knots in the BL Lac objects have flux within 
25% to 100% of the flux of the core, significantly brighter than fast knots in these objects. 

If fast knots represent forward shocks and slow knots correspond to reverse shocks, 
then the forward shock is stronger in the quasars and, perhaps, the radio galaxies, while the 
reverse shock dominates in the jets of BL Lac objects. A reverse shock will be strong relative 
to its corresponding forward shock when the disturbance is prolonged such that the faster 
flow enters the rear of the shock structure over an extended period of time. Otherwise, only 
the forward shock will be prominent and the knot will have a Lorentz factor only slightly 
(<6%) less than that of the forward shock front (e.g., Sokolov, Marscher, & McHardy 
2004). BL Lac objects might therefore have more prolonged disturbances of lower amplitude 
than those in quasars. If jets possess spine-sheath structure (Laing et al. 1999), then an 
alternative explanation for the above effect could be that the power of high-r spines in 
BL Lacs is lower than in quasars. The proposal that the nature of shock waves in quasars 
is different from that of BL Lac objects has been suggested by Wardle et al. (1994), who 
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interpreted the evolution of the flux and polarization of the quasar 3C 345 as forward shocks 
that have greater speed than the underlying jet. This is in contrast to the BL Lac object 
OJ 287, studied by Cawthorne & Wardle (1988), who modeled the polarization, kinematics, 
and X-ray variability as reverse shocks, which are slower than the underlying jet. Our data 
support the idea that this could be a key difference between the two classes of blazars. 

6.6. Trailing Shocks 

Numerical hydrodynamical simulations (Agudo et al. 2001; Aloy et al. 2003) show that 
the interaction of disturbances in the jet flow with the underlying jet and/or the external 
medium can play a significant role in the variability of the jet emission. In particular, multiple 
conical shocks can form behind a strong shock wave propagating down the jet. These trailing 
shocks appear to be released in the wake of the primary superluminal component rather 
than ejected from the core. The simulations show that the ratio between the apparent 
velocity of the main and trailing components is a function of distance along the jet: the 
closer to the core a trailing component first appears the slower it moves, although for a given 
trailing component the simulations predict deceleration with distance from the core. Trailing 
components are oblique shocks in three-dimensional models and should possess different 
polarization properties than the leading shock. Gomez et al. (2001) identified several knots 
in the jet of the radio galaxy 3C 120 as having the characteristics of trailing shocks. Our 
sequences of high resolution VLBA images reveal a number of jet features in different sources 
that have properties matching those expected for this phenomenon. These are (see Fig. 16): 
cl, c2 behind knot CI in the radio galaxy 3C 111; o2 behind ol found previously by Gomez et 
al. (2001) in 3C 120; 61, 62 following Bl, B2, respectively, in the quasar 3C 273; c9 behind C9 
in 3C 345; 65, 66 following B5, B6, respectively, in CTA 102; and 63 behind B3 in 3C 454.3. 
Components c2 (3C 111), 62 (3C 273), and 65 and 66 (CTA 102) emerge from bright knots 
at some distance from the core (see Fig. 16), while the remainder are already trailing bright 
features from the first epoch of our observations. All of them have apparent speeds less 
than the corresponding main component and exhibit different EVPAs than the fast feature 
(see Fig. 17). Although our observations do not allow us to check every trailing component 
for deceleration, a decrease in apparent speed is pronounced for cl and o2 in 3C 111 and 
3C 120, respectively, and can be inferred for 65 in CTA 102 (see Fig. 33). For 65 in CTA 102 
a second-order polynomial does not satisfy the x 2 test but improves goodness-of-fit to the 
data significantly. 

We characterize each trailing component by three parameters: /3 app is the average appar- 
ent speed listed in Table 5, P&pp/Pmp is the ratio of apparent velocities of trailing component 
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to corresponding main component, and R is the average deprojected distance from the core 
at which a trailing component is detected in our data (R is calculated using < R > indicated 
in Table 5 and jet parameters given in Table 11). Figure 34 (left panel) demonstrates that 
trailing components form in the wake of bright knots having values of T extending from 5 to 
>20 and different values of (3 app /f3^ (all of them <1) are observed for similar Lorentz fac- 
tors. Figure 34 (right panel) shows that there is an increase of the apparent speed of trailing 
components with distance from the core (coefficient of correlation 0.71). According to Agudo 
et al. (2001), trailing components represent pinch waves excited by the main disturbance, 
and an increase of their speed at larger distance reflects acceleration of the expanding jet. 
Our data shown in Figure 34 imply such an acceleration of the underlying jet that appears 
not to depend strongly on the Lorentz factor of the main disturbance, consistent with the 
conversion of internal energy into bulk kinetic energy that accompanies expansion. 

Figure 34 (left panel) shows a possible inverse correlation between [3 app / P^, p and Lorentz 
factor of the main disturbance (coefficient of correlation —0.40). The correlation is most likely 
the result of selection effects. Detection of a trailing component close to a leading shock with a 
high Lorentz factor is complicated by stretching of the longitudual size of the main component 
in the observer's frame by a factor oc T owing to light-travel time delays. Detection of trailing 
components that lag greatly behind a main disturbance with a low Lorentz factor is hampered 
by difficulties in associating with confidence trailing components to the leading knot. For 
this reason, many of the subluminal or quasi-stationary features detected near the core (A- 
components) in the majority of sources in the sample could represent a superposition of many 
trailing components formed behind a number of superluminal features. This possibility is 
valid only if the direction of the jet "nozzle" changes between ejections, since a second major 
disturbance passing through a trailing shock would destroy it. 

Kellermann et al. (2004) have found that there is a systematic decrease in (3 app with 
increasing wavelength, which they suggest results from sampling different parts of the jet 
structure at different frequencies. Although this might be the case for average apparent 
speeds derived from surveys at different wavelengths, for individual sources with superluminal 
components detected at similar distances at different frequencies, the effect can be caused by 
the structure containing the leading perturbation plus slower trailing shocks being unresolved 
at the longer wavelengths. Perhaps the difference in the proper motions of the polarized and 
total intensity components B8p and B8 noted in the quasar 0528+134 (see Table 5) is an 
example of this effect, with the polarized intensity image (which reveals the leading compact 
component) playing the role of the finer resolution of shorter wavelength observations. 
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6.7. Frequency of Superluminal Ejections 

Our monitoring is unique in terms of the number of blazars observed in a regular manner 
at high angular resolution over 3 years. The excellent time coverage allows us to determine 
the rate of superluminal ejections, thought to be controlled by activity in the cental engine. 
We list in Table 11 the rate, / e j (multiplied by the time dilation factor 1+z), of superluminal 
ejections over three successive years based on the results given in Table 5. In Figure 35 
we plot the ejection rate versus the average jet Lorentz factor for blazars. There is a trend 
suggesting a positive correlation between / e j and T (coefficient of correlation 0.5) The formal 
t^ v test rejects the hypothesis that there is no relation between / e j and T at a confidence level 
£=0.05. However, the trend could be an artifact of the bias toward highly Doppler boosted 
objects in our sample. Nevertheless, blazars form the largest class of identified 7-ray sources 
(Hartman et al. 1999), for which a connection between 7-ray events and radio jet activity 
has been found (Jorstad et al. 2001b). If nonthermal flares are associated with ejections 
of new superluminal knots, as seems to be the case at mm wavelengths (Savolainen et al. 
2002), then we would expect 7-ray light curves to reflect the rate of ejections. In this case, 
the possible correlation between the Lorentz factor and ejection rate would imply that the 
most variable 7-ray sources (in terms of major flares per year) should possess the highest 
Lorentz factors. This can be tested by 7-ray light curves of blazars obtained by the Gamma 
Ray Large Area Space Telescope (GLAST, expected to begin operation in 2007), along with 
VLBI monitoring of the radio jets. 



7. Summary 

In this paper we have presented the entire data set collected during a three-year program 
of monitoring AGNs using the VLBA, and have discussed the main results relating to the 
kinematics of jets. The sequences of images reveal short timescales (for some sources shorter 
than 2 months) of variability of the jet structure and even more rapid variability in polar- 
ization. This program illustrates the importance of intensive monitoring for understanding 
jet physics in superluminal radio sources. 

We have measured the apparent speed of 106 features in the inner jets (within 4 mas 
of the core) of two radio galaxies, five BL Lac objects, and eight quasars from 1998 March 
to 2001 April. Superluminal apparent speeds occur in 80% of the knots, 26% of which 
show statistically significant deviations from ballistic motion. The majority of non-ballistic 
components undergo an increase of apparent speed with distance from the core, although 
local decelerations are observed in some cases alongside the general acceleration. This could 
be the result of physical accelerations or from selection of sources whose angles to the line of 
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sight < sin _1 (l/r) near the core and closer to this value farther out. Many of the jets contain 
both very fast and much slower moving features, which might be explained as forward and 
reverse shocks, respectively. It appears that fast features are pronounced in the jets of the 
quasars, while slower moving knots dominate in the BL Lac jets. This suggests a different 
nature of the main disturbances seen in the jets of the two classes of blazars. 

The properties of 11% of the superluminal components are consistent with the charac- 
teristics of trailing shocks expected to form in the wake of strong disturbances in the flow. 
Trailing components emerging from bright knots farther down the jet have a faster apparent 
speed than those generated near the core while a given trailing feature most likely decelerates 
along the jet. According to the numerical simulations by Agudo et al. (2001), the trailing 
components are caused by the triggering of pinch modes by the main disturbance. The jet 
accelerates close to the core as internal energy is converted into bulk kinetic energy, leading 
to trailing shocks being faster farther downstream. 

In four sources (3C 120, OJ 287, 3C 345, and CTA 102) the behavior of jet components 
appears to be affected by the interaction with the external medium (see §A). For each source 
a number of components decelerate, change trajectory, brighten in total and polarized flux, 
and undergo a rotation of the EVPA at the same distance from the core along an edge of 
the jet. This implies the existence of gas clouds a few parsecs to >1 kpc (deprojected) from 
the central engine, intermediate between the locations of the very dense clouds of the broad 
line regions and more rarefied clouds of the narrow line regions. 

Using measurements from the VLBI images of parameters such as flux density, apparent 
speed, and size of components, we have estimated Lorentz and Doppler factors as well as 
viewing angles for superluminal knots and the opening angle of each jet. This is a new 
method to define jet parameters, based on the assumption that the decay in flux of the 
superluminal components is caused by radiative losses rather than by cooling from expansion, 
and is subject to light-travel delays. We demonstrate that at high radio frequencies these 
assumptions are most likely correct. The derived parameters of the jets indicate that in 
our sample the quasars have the highest Doppler factors (5) and smallest viewing (6 Q ) and 
opening (6) angles, while the two radio galaxies possess significantly lower Doppler factors, 
larger viewing angles, and wider opening angles despite their "blazar-like" radio properties. 
This implies that in the 3-dimensional parameter space (S,Q o ,0) the radio galaxies, BL Lacs, 
and quasars in our sample occupy different regions. The regions of the quasars and BL Lac 
objects partly overlap, while the radio galaxies are significantly distinguished from the blazars 
in all three parameters. Since our sample is not a complete one, the major differences 
suggested by these results need to be investigated further with larger, complete samples. 
The inferred relationship between the half opening angle and the Lorentz factor agrees with 
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the expectation of gas-dynamic models that predict smaller values of 9 for higher Lorentz 
factors and a dependence of 9 on the ratio of the external and internal jet pressure. The 
best approximation to the relation is very close to P cxt /P ~ 1/3. This does not, however, 
exclude collimation by magnetic pinching, which might produce a similar 9 — T relation. 

We have estimated the intrinsic brightness temperatures of jet components in the 
quasars, BL Lacs, and radio galaxies on parsec scales, obtaining averages of 1.1 x 10 9 K, 
5.5 x 10 7 K, and 3.5 x 10 9 K, respectively. Comparison of these values with the equipartition 
brightness temperature of the optically thick part of the jets, Tb 5 i nt = 2 — 5 x 10 10 K, suggests 
a stronger magnetic field in the BL Lac objects. 

There is a possible positive correlation between the Lorentz factor of the jet and ejection 
rate of superluminal components for the blazars. This can be tested with the 7-ray light 
curves that will be measured by the GLAST mission. 

A. Notes on Individual Sources 

3C 66A: The radio jet consists of a bright, polarized core with weak, diffuse extended 
structure (Fig. 1). The set of images allows us to trace components CI, C2, C3, and C4 
at many epochs (Fig. 16). They are characterized by moderate velocity, 2-5c, typical of 
those measured in BL Lac objects at lower frequencies (Gabuzda et al. 1994). Component 
CI is probably feature C from Jorstad et al. (2001a), where it was thought to be stationary. 
The centroids of Al and A2 appear to move toward the core at mildly superluminal speeds; 
however, these are extended features, hence the shifts in position could be due to changes 
in their brightness distributions. This is supported in the case of A2, whose proper motion 
makes an angle of 45° to the line between the core and the feature's centroid. Figure 16 shows 
possible identifications of very weak components BA, B6, and B7. The designation of knots 
follows that of Jorstad et al. (2001a); Figure 16 includes their data, marked by triangles. B 
knots move at a high apparent speed, > 20c. The reality of this faster motion is supported 
by a region of enhanced polarized intensity that progresses rapidly downstream at the same 
epochs as when the high apparent speed is detected in the total intensity images (see Fig. 
36). It appears that the more slowly moving C components and highly polarized, fast B 
components coexist in the jet. One of possible physical explanations for this dichotomy is 
a double shock structure, with the faster forward shock being weaker than slower reverse 
shock (see §6.5). 

3C 111: The innermost jet of this galaxy exhibits strong activity. At least 7 moving 
components (Bl — B7, see Fig. 16) can be identified inside 1 mas of the core. These knots 
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move ballistically at a typical apparent speed ~ 3c. There are two bright features near the 
core, Al and A2, that are stationary according to the x 2 test. The most prominent feature 
of the jet is an extended, polarized blob northeast of the core. The feature appears to be 
associated with a major outburst observed at 90 GHz in 1996 January (Alef et al. 1998), 
at 37 GHz in 1996 March, and at 22 GHz in 1996 August (Terasranta et al. 2004). The 
first ejection around this time appears to have occurred earlier than the derived epoch of 
ejection T (see Table 5), since complex structure of the core region was detected already in 
a 43 GHz VLBI map obtained in 1996 July (Alef et al. 1998). This implies slower motion of 
the disturbance near features Al and A2, which agrees with the description of the structure 
by Alef et al. (1998) as "a mix of dimming and flaring of stationary and moving emission 
components." In our observations we identify the brightest part of the outburst remnant 
as CI, which moves ballistically at a higher apparent speed than that of the innermost 
components. Knots cl and c2 detected behind Cl seem to belong to the same disturbance. 
At the last 3 epochs the parameters of Cl, cl, and c2 are obtained after tapering the images 
by reducing the weighting of the data from the longest baselines, and are not used to compute 
the jet velocities. 

0420—014: The parsec-scale jet of this quasar is strongly dominated by the VLBI core, 
which has a complex polarization pattern (see Figs. 3 and 17), implying multi-component 
substructure in the core region. We detect two moving components in the innermost part of 
the jet. Knot Bl has slower motion near the core, ~7c (Fig. 30), and accelerates significantly 
beyond 0.2 mas as its trajectory turns from —100° to —175° (Fig. 18). The bright component 
B2 traced inside 0.15 mas from the core moves ballistically along ~ —72° at a similar 
apparent speed as that derived for Bl near the core. The behavior of the knots is consistent 
with the common curved trajectory for jet components found by Britzen et al. (2000). 

3C 120: Our observations continue to trace components d, h, I, ol, and o detected 
in Gomez et al. (2001), and we identify new components t, u, ul, and v (Fig. 16). Com- 
ponents d and h move ballistically with the same apparent speed found previously, while 
/ has decelerated since earlier epochs. Components ol and o2 represent the front section 
of the major disturbance designated in Gomez et al. (2001) as complex o. Component ol 
brightens in total and polarized flux at ~3.4 mas, similar to the behavior seen for h and I 
(Gomez et al. 2001). After flaring, ol and o2 diverge in their behavior: ol accelerates and 
o2 decelerates (Fig. 16). Such a development is expected if o2 is a trailing component that 
forms in the wake of the accelerating knot ol, which represents the major disturbance in the 
flow (Gomez et al. 2001). The most prominent feature at later epochs is t, which evolves 
in the same manner as I. The components have nearly equal proper motions, 2.17±0.06 (/, 
Gomez et al. 2001) and 2.04±0.08 mas yr _1 (t), and similar trajectories (Fig. 18). They 
both undergo a flare in total and polarized flux at a location between 2.5 and 3.5 mas from 
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the core and the EVPAs of both rotate during the flare (Fig. 37). Gomez et al. (2000) 
explained these properties as the results of interaction between the jet and an interstellar 
cloud. The behavior of component t is consistent with this, which suggests that it interacts 
with the same cloud. However, there is a significant shift (~ 0.6 mas) between the projected 
positions of the peaks of the flares (see Fig. 37). Most likely, the shift is a consequence of 
the slightly more southern trajectory of t (Fig. 18). For a viewing angle ~ 21° and opening 
angle ~ 4° (see Table 11), the shift by 0.6 mas gives a lower limit to the size of the cloud of 
1 pc at a deprojected distance from the VLBI core of ~7-8 pc. 

0528+134: Components B3 and B4 detected in Jorstad et al. (2001a) are seen during 
the early epochs of our monitoring and four new components are ejected later (Fig. 16). 
All components have similar trajectories, with a sharp bend to the north beyond 0.5 mas 
from the core, as previously noted by other authors (e.g., Britzen et al. 1999). The projected 
trajectories of the knots are shifted relative to each other (see Fig. 18), which suggests 
different angles of ejection. This is supported by a difference in the apparent speed near the 
core, from > 10c for the fastest and most eastern components to ~ 3c for the most northern 
component. We find progressive acceleration in at least three components (Fig. 30). The 
most northern knot, B8, which is strongly polarized and the brightest feature in the jet (Fig. 
17), shows significantly different motion of the total and polarized intensity peaks. The 
polarized region, B8p, moves ballistically ahead of the main concentration of total intensity, 
B8. Table 5 lists the parameters of B8 and B8p separately. We suggest that B8p represents 
a leading shock of a strong disturbance with complex structure (see §6.6). 

OJ 287: The jet appearance changes considerably during the observations (Fig. 6), 
largely due to brightness variability of three prominent jet features, Al, A2, and A3 (Fig. 17). 
On average, Al and A2 separate from the core at subluminal speeds, although significant 
fluctuations relative to linear motion are observed (see Fig. 38). Component A2 might be 
associated with CI, classified by Jorstad et al. (2001a) as a stationary feature. Table 5 gives 
the parameters of A2 based on our data alone, and also based on these data combined with 
those of Jorstad et al. (2001a) under the assumption that A2 and CI are the same feature. 
The latter is listed as A2* in Table 5. The addition of the previous data decreases the 
proper motion by a factor of ~2 and in both cases the proper motion does not exceed zero 
at the 2a level. This suggests that the feature may fluctuate about a stationary position 
rather than translate downstream. We identify several fast moving components in the jet. 
Knots BA — B7 (the designation continues the scheme of Jorstad et al. 2001a) have very 
short lifetimes of ~1 yr at 43 GHz, and can be resolved only after they have passed A2 (see 
Fig. 17). However, a notable brightening of A2 occurs at the epochs when a B component 
likely reached A2 (see Fig. 16). The B components show high superluminal speeds that are 
slightly slower than those of knot K3/U3 in Homan et al. (2001). Knot K3/U3 is observed 
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to accelerate and its motion at 22/15 GHz is measured ~2 times farther from the core than 
the locations of the B components in our images. 

A distinguishing feature of the jet is component A3, the formal solution for the motion 
of which is upstream, toward the core. This motion is especially pronounced in the likely 
event that A3 is associated with C2 identified by Jorstad et al. (2001a), in which case A3 
also undergoes a significant change in position angle relative to the core (Fig. 18). The 
latter can be explained by swinging/precession of the jet nozzle or by a translation of the 
core, either of which agrees with changes in trajectory of components near the core from 
~ —90° for Bl (Jorstad et al. 2001a) to ~ —115° for B6. The alignment between the 
innermost jet direction and the position of component A3 suggests that the feature is the 
result of interaction between the jet flow and the external medium. The upstream apparent 
motion of A3 might then trace the location of an interstellar cloud (projected on the sky) 
when the intrinsic position angle of the jet changes. In this case, the dramatic fading of the 
feature toward the end of our observations (see Fig. 6) could be connected with an intrinsic 
change of the inner jet direction (component B6 has the most southern projected position 
angle during 1995-2000) so that the jet flow misses the cloud. This interpretation agrees 
with the findings of Tateyama & Kingham (2004). These authors used the Radio Reference 
Frame Image Database data and determined that the projected position angle of the inner 
jet of OJ 287 rotated clockwise by ~ 30° in 8 yr, in agreement with a ballistic precession 
model if the velocities of the components decrease beyond ~1 mas from the core. We suggest 
that the interaction with an extragalactic cloud might cause the deceleration of components. 
From the derived parameters of the jet (see Table 11), we place the cloud at a deprojected 
distance of ~ 90 pc from the core. Acceleration of K3/U3 detected by Homan et al. (2001) 
takes place between 0.2 and 1 mas from the core and reflects an increase in the apparent 
speed beyond Al and A2 and before A3 located at ~1 mas. 

3C 273: The most prominent features in the images (see Fig. 7) are an unpolarized core 
and strongly polarized moving knots Bl and B2 that both brighten at ~0.8 mas (Fig. 39). 
Each component becomes strongly polarized just as it leaves the core region (~0.3 mas from 
the core, Fig. 17) while the core has a low polarization (comparable with the noise level). 
The similarity of the EVPA and fractional polarization at 7 (brighest polarized feature), 3, 
and 1.3/0.85 mm suggests a low fractional polarization in the core at the high frequencies, 
where Faraday effects are not important. This implies that the core has intrinsically low 
polarization, which might be caused by a strongly turbulent magnetic field on scales smaller 
than the 7 mm synthesized beam. After they flare, Bl and B2 fade dramatically and expand 
both along and transverse to the jet direction. This results in the formation of new features 
bl and b s related to SI, and 62 related to B2 (Fig. 16). The b components might form in 
the interaction of relativistic shocks connected with features Bl and B2 and the underlying 
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jet flow. These components should have lower apparent speed than the main disturbances 
(Agudo et al. 2001), but this is not the case for b s . Alternatively, the B components are 
fairly extended features in a broad jet and the b knots might represent sections of the jet 
where the velocity vectors lie at different angles to the line of sight. The latter interpretation 
is supported by a difference in the projected trajectories of the B and b components (Fig. 
18) and the strong dependence between the projected position angle and apparent speed of 
components, with the more northern knots being slower (Fig. 40). A substantial velocity 
gradient across the jet is also expected according to a model where the underlying jet flow has 
a double-helix structure and a slower Lorentz factor than that of disturbances (Lobanov & 
Zensus 2001). The brightening of Bl and B2 at ~0.8 mas from the core might be connected 
with either (1) a change of the jet direction (a slight curvature of the projected trajectories 
of both components is seen in Fig. 18 at this distance), (2) an intensity peak of a threadlike 
pattern of the underlying flow where the local velocity vector bents toward the line of sight 
(Lobanov & Zensus 2001), or (3) interaction with the external medium. However, a stable 
polarization direction in B2 before, during, and after the flare (see epochs 1999.76, 1999.93, 
and 2000.07 in Fig. 7) does not support (3). 

At many epochs the images contain a diffuse, highly polarized (up to 50%) feature, CI, 
farther down the jet (see Fig. 41), moving at an apparent speed of ~ 7c. According to the 
time of ejection, CI corresponds to knot Gl found by Jorstad et al. (2001a). If CI and Gl 
are the same feature, it has maintained a nearly constant proper motion over ~5 yr. 

3C 279: Our data reveal motion of at least 8 components, C9 — C16 (our designations 
continue those of Wehrle et al. 2001), with apparent speeds up to 17c, three times higher 
than previously reported velocities measured since 1991 (Wehrle et al. 2001; Homan et al. 
2001; Jorstad et al. 2001a). The identification of components in Figure 16 includes results at 
43 GHz from Wehrle et al. (2001) (triangles) and Jorstad et al. (2001a) (squares). Compo- 
nents C9 — C16 are shown in Figure 17 at the last epoch of monitoring. Jorstad et al. (2004) 
suggest that the remarkable change in apparent speed is caused by a shift in the direction of 
the jet "nozzle" while the Lorentz factor remains constant. This conclusion is supported by 
a gradual decrease of the projected position angle O of each successively ejected component 
(starting with C9). In this interpretation, the jet flow is characterized by T > 20 while O 
changes from 0.5° for components with the lowest speed (5.5c for C8 found by Wehrle et al. 
2001) to 1.6° for components with the highest speed. It is also possible that component C8 
accelerates (see Fig. 16) due to an increase in the Lorentz factor of the flow. We designate 
this possible component as C8*; it corresponds to (78 before the acceleration and to C9 after 
the acceleration. As indicated in Table 6, the change in the velocity component parallel to 
jet, is significant. Another method, which we apply in this paper to estimate Doppler 
and Lorentz factors and viewing angle for each superluminal knot (see §6), suggests that 
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both the Lorentz factor and the angle between the jet axis and line of sight for C9 — C16 are 
different from those of C8. In any event, these results are generally consistent with viewing 
angle playing the major role in the change of apparent speed. 

Contemporaneous to these changes in the inner jet, an increase in the apparent speed 
is observed for CA at ~3 mas from the core at epoch ~1998.2 (Homan et al. 2003). Our 
data, combined with previously published results at 43 GHz (Wehrle et al. 2001; Jorstad et 
al. 2001a), show that variations in the trajectory and apparent speed of CA occurred at least 
twice between 1995 and 2001. This identifies three segments (a, b, and c) of the trajectory 
(denoted by dotted lines in Fig. 18) where the motion is ballistic according to a x 2 test. The 
ballistic sections of the trajectory define two epochs, 1996.1±0.1 and 1998.07±0.03, when the 
apparent speed changes. The values of /3 app corresponding to these paths are given in Table 
5 using designations CAa, CAb, and CAc, respectively. Table 5 shows that CA had a high 
apparent speed ((3 apP ~ 13c) in 1995 and decelerated dramatically in 1996, then accelerated 
again from 1998 until at least 2002 according to Homan et al. (2003). Both periods of 
acceleration correspond to fading of the component, while during times of lower apparent 
speed CA brightened. Such behavior may be the result of helical motion, with the radius of 
the helix increasing outward. 

1510—089: The radio emission of this quasar is strongly core dominated with usually 
weak extended jet structure. Our data indicate that the quasar possesses an ultra-relativistic 
jet with apparent speed up to ~46c. However, even bimonthly observations are not sufficient 
to state that the identification of components shown in Figures 9, 16, and 42 is unique. On 
the other hand, this identification is supported by the total and polarized intensity images 
along with the total intensity light curves at 37 GHz (Terasranta et al. 2004) and 14.5 GHz 
(UMRAO data base). We detect three moving components, Bl (Fig. 9), B2, and B3 (Fig. 
42). In all three cases we observe a similar polarization behavior in the core region when 
the component is emerging from the core - the core region is strongly polarized with EVPA 
along the jet direction (see Fig. 9). In the cases of components Bl and -B3 polarized features 
with similar EVPA are subsequently observed propogating down the jet (Fig. 42). For all 
three components the time of the ejection falls on the rising branch of an outburst in the 
total intensity light curves (an especially prominent flare is connected with the ejection of 
component B2). The modeling of the jet structure shows two features, Al and A2, that 
fluctuate near positions 0.14 and 0.54 mas, respectively (Fig. 16). Component Al is usually 
canot be resolved from a moving component at epochs when the latter is close to Al, while 
A2 shifts downstream when a moving component is approching and returns back when the 
moving knot has passed by (see Fig. 16). Such a behavior of jet features formed by stationary 
shocks is seen in numerical simulations of jets with variable input flow velocities (Aloy et al. 
2003). At many epochs, a feature, S, with flux significantly above the noise is seen southeast 
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of the core (Fig. 17), close to the position angle of the arcsecond scale jet (Homan et al. 
2002b). On the other hand, this component is located at a similar distance from the core 
as Al but in the opposite direction. This suggests possible artificial generation of the knot 
owing to modeling of an actually slightly extended core using a point-like source. However, 
there is no correlation between the positions and fluxes of S and Al. The position angle of 
S varies within the range 71° to 168° and the flux varies from 0.03 Jy to 0.5 Jy, while the 
values of 6 and the flux of Al are fairly stable (parameters of Al and A2 in Table 5 are 
computed for epochs when the features are not confused with moving knots). Differences in 
the properties of S and Al suggest that these features along with AO constitute a complex 
core region inside the 43 GHz synthesized beam. 

3C 345: At many epochs the core of this quasar is not the brightest feature of the 
jet (Fig. 10). The core is usually unpolarized at 43 GHz (Fig. 17), although the fractional 
polarization of the component nearest to the core, Al, can be as high as 20%. The latter fact 
suggests the low polarization of the core is intrinsic rather than caused by Faraday effects. 
We detect seven superluminal knots with average apparent speeds in the range 9-24c (Fig. 16, 
where the component designations follow Ros et al. 2000), and at least four of these exhibit 
acceleration or deceleration (Figs. 18 and 30). For knots C10 — CT3 statistically significant 
acceleration is found only for Cll. Similar acceleration was reported for components C4 
and C7 by Biretta et al. (1983) and Ros et al. (2000), although beyond 2 mas from the core 
C7 decelerates as indicated by both (x\\ and interpolation of its position relative to 1996.81 
as observed by Ros et al. (2000) (asterisks in Fig. 16). It appears that C8 repeats the stages 
of C7, although it experiences an earlier and greater deceleration, possibly signifying a more 
complex trajectory (see Fig. 18). The change in apparent speed and the twisted trajectory 
are also pronounced for C9. Although the historically extrapolated position of C9 (using the 
average apparent speed of 15.6c) agrees well with the measured position at epoch 1996.81 
(Ros et al. 2000), it decelerates significantly from ~ 20c near the core to ~ 10c at 1 mas, 
and then accelerates beyond 1.5 mas. The deceleration is accompanied by brightening of the 
knot and rotation of the EVPA by ~ 50°, making it almost perpendicular to the jet direction 
(see epochs 1998.76 and 1998.94 in Fig. 10). In addition C9 has the most southern position 
angle in the jet during the deceleration period (6 ~ —116°), similar to that observed for C8 
by Ros et al. (2000, 6 ~ -122°). Wardle et al. (1994) interpreted a change they observed 
in apparent speed in terms of a flow with constant Lorentz factor along a bent trajectory. 
Component C9 expands as it propagates down the jet and leaves behind a complex structure 
that consists of knot c9 moving at apparent speed ~ 9c (Fig. 16) plus two extended features, 
c9n and c9s, that form later at the northern and southern edges of the jet, respectively (Fig. 
17). The majority of moving features have a trajectory with zig-zags in the north-south 
direction (although the net motion is to the west), an increase in flux, and rotation of the 
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EVPA at a distance ~0.8-l mas from the core (Fig. 43). These variations could be the result 
of interaction with the external medium, as suggested by Ros et al. (2000). We estimate the 
deprojected distance of this possible interaction as ~300 pc. 

1803+784: The curved jet of this BL Lac object has been interpreted as the result of 
helical structure (Britzen et al. 2001; Tateyama et al. 2002). Our high resolution 43 GHz 
images show at least two superluminal components near the core, Bl and B2. These follow 
a similar curved path (Fig. 18) and have average apparent speeds of ~ 16c, with acceleration 
and deceleration along the trajectory (Fig. 30). A passage from deceleration to acceleration 
and vice versa is detected three times for Bl, and locations of these transitions along the jet 
coincide with bends in the trajectory (compare Figs. 18 and 30). This behavior supports 
the hypothesis of helical structure of the jet. Gabuzda & Chernetskii (2003) classify the jet 
feature ~1.5 mas from the core as a stationary knot, S. In our data this feature appears to 
be associated with A2, and has an upper limit to its apparent speed of ~ 1.5c (Fig. 16). 
No motion is detected for bright knot Al near the core. Tateyama et al. (2002) found that 
a superluminal component at 8.5 GHz only becomes visible at a distance >1 mas from the 
core. Moving component CI is detected beyond Al and has an apparent speed of ~ 10c, 
similar to that found by Britzen & Krichbaum (1995). The fractional polarization in the jet is 
high. For example, A2 is polarized up to 30% with the EVPA aligned with the jet direction. 
In contrast, the EVPA in the core region is essentiallyradial (see Fig. 17), supporting the 
inference of a toroidal magnetic field geometry Gabuzda (2000), but only if the viewing angle 
is almost 0° in the immediate vicinity of the core. It is interesting to note that the viewing 
angle of the jet derived in §6 is indeed quite small. 

1823+568: The 43 GHz images of this BL Lac object contain bright stable features 
A1,A2, and A3 within 0.5 mas of the core (Fig. 16). All of these components are strongly 
polarized, especially A3, where the fractional polarization reaches 25% (Fig. 17). Four 
components, A2, A3, C2, and CI, detected in our images can be identified with components 
K4, K3, K2, and Kl, respectively, from Gabuzda & Cawthorne (2000). The locations of 
A2 and A3 are slightly farther from the core than K4 and K3 (epoch 1995.58). This can 
be attributed to a shift in the core position at 43 GHz relative to 22 GHz due to opacity 
effects, in which case the positions of these features would be stable over a period of 5 years. 
Despite excellent time coverage (Fig. 12), we do not detect motion inside 1 mas of the core 
in our observations. This has several possible explanations: 1) ejections of superluminal 
components are rare in this object, with mean time between ejections > 20 years in the 
observer's frame (see Table 5); 2) moving components are much fainter than stationary 
features; or 3) the innermost part of the jet is very close to the line of sight. However, we 
do detect two moving components, CI and C2, beyond 1 mas from the core where the jet 
bends slightly to the west (see Fig. 17). Knot CI has a low apparent speed of 3.8c that 
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appears to be essentially constant, since the linear extrapolation of its motion back to epoch 
1995.58 agrees well (within 0.03 mas) with the position of knot Kl in the image of Gabuzda 
& Cawthorne (2000). For the fastest moving feature in the jet, C2 (f3 app ~ 8c), extrapolation 
from our sequence produces a position of the component at epoch 1995.58 that is farther 
from the core by ~0.1 mas than was observed for K2. This implies a slower apparent speed 
farther from the core. Motion of C2 is not detected in our images after it passes CI (Fig. 
16). 

BL Lac: Our data for BL Lac have been analyzed in Stirling et al. (2003). At least 
three new components are ejected during the monitoring program (Fig. 16). These, together 
with previously ejected knots observed by Denn, Mutel, & Marscher (2000), form a sequence 
of superluminal features that appear to emerge from the core ballistically along different 
position angles, indicating a change in direction of the jet "nozzle." The trajectories of 
the knots exhibit curvature that is especially pronounced for components 510, 511, and 
513 (Fig. 18). Despite different average apparent speeds, the components accelerate in a 
similar manner suggestive of a small change in the Lorentz factor and a substantial change 
in direction (see Table 6 and Fig. 30). The core region can be represented by two point-like 
features, CI and C2 (Stirling et al. 2003, and Fig. 17). In the inner jet two bright stationary 
knots, Al and A2, are observed at all epochs (see Figs. 16 and 17). 

CTA 102: The apparent velocity of jet components in this quasar is fairly high, j3 app ~ 
18c, in agreement with previous findings (Jorstad et al. 2001a; Rantakyro et al. 2003). We 
label newly detected moving components as B5 and B6 (Fig. 16) following the scheme of 
Jorstad et al. (2001a) and rename J in Rantakyro et al. (2003) as B4 based on its estimated 
ejection time. Components B5 and B6 have ballistic motion. However, at ~ 0.2 mas from 
the core both knots split, forming more slowly moving components 65 and 66 behind the 
main features. At some epochs (e.g., Fig. 17) Bb is significantly brighter than the core both 
in total and polarized intensity. There are two features, Al (near the core) and C (~2 mas 
from the core), that do not move according to a x 2 test. Component C was detected by 
Jorstad et al. (2001a) and classified as probably stationary, although Rantakyro et al. (2003) 
associated the feature with different moving components. As noted by Jorstad et al. (2001a), 
C is strongly variable in brightness and located where the jet bends sharply to the south. In 
1998-1999, C was bright, very diffuse, and showed a double structure that includes knot c 
closer to the core (see Fig. 17). We detect no motion of c (Fig. 16). Both components fade 
dramatically by the end of the monitoring (Fig. 14) despite the approach of the growing inner 
part of the jet toward C. This behavior appears similar to that observed for component A3 
in OJ 287 (see above), where we associate the fading of A3 with a change in the jet direction. 
Figure 18 shows the trajectories of B5 and B6, as well as Bl and -B3 from Jorstad et al. 
(2001a), plus the position of C during both sets of observations. The location of C fluctuates 
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over the range 1-1.5 mas from the core in right ascension and within 1.4-1.6 mas (to the south) 
in declination. The innermost jet shows a significant (~ 20°) difference in projected position 
angle of B5 and B6 relative to Bl — B3, suggestive of wobbling of the jet "nozzle" and/or 
change of the VLBI core position. These properties imply interaction between the moving 
flow and external medium at a deprojected distance ~1.4 kpc from the core that results in 
the formation of a diffuse feature. 

3C 454.3: This quasar has very complex structure at 43 GHz. Figure 16 plots only 
bright components (flux density >0.1 Jy). One of the most interesting features of the jet 
is highly polarized component C. At all epochs (see Figs. 15 and 17) the EVPA of the 
brightness peak of C aligns with the jet axis. No change in separation of the component 
from the core has been seen during either previous observations (Jorstad et al. 2001a) or our 
current sequence of images. However, combining the data from the two studies reveals a drift 
of C to the southwest (see Fig. 18) at an apparent speed of ~ 1.6c between 1995 and 2001 
(the component defined in this way is referred to as C* in Table 5), while the size of the knot 
remains fairly stable at ~ 0.25 ± 0.10 mas. The absence of proper motion within each data 
set suggests that the observed difference in position of the feature (~0.2 mas in declination 
and ~0.1 mas in right ascension) could be caused by a displacement to the northeast of 
the VLBI core between 1997.6 and 1998.2. Beyond C, trajectories of moving components 
separate into two groups: northwestern (B3,B6) and southwestern (-B4), with a difference 
between the projected trajectories that can reach ~ 70° (Fig. 18). This is consistent with 
the direction of the jet near the location of C lying very close to the line of sight, thereby 
producing large differences in projected trajectories of moving components beyond C even if 
the intrinsic trajectories differ only slightly. The proper motions that we measure cover the 
entire range from 0.12 to 0.53 mas yr _1 found previously in the jet of 3C 454.3 (Gomez et 
al. 1999; Jorstad et al. 2001a, and references therein). Component B6 accelerates markedly 
(see Table 6). A similar acceleration is derived for a bright feature detected in 1998 beyond 
C if we associate it with B3 found by Jorstad et al. (2001a) based on its brightness, position, 
and time of ejection. This gives an increase in apparent speed of B3 from 6c near the core 
to 13c beyond C. Both of these components belong to the northwestern group (Fig. 18). 
Northwestern components appear to accelerate significantly without a change in direction of 
motion while southwestern knots (e.g., BA) show mostly a change in direction (Fig 30). A 
diffuse feature (designated in Fig. 16 by the symbol +) is located in the northern part of the 
jet and appears to move after B4 passes C. Figure 17 displays component S to the south 
of C when B4 approaches C. As observed in other objects discussed above, more slowly 
moving knots 63 and b6 trail some of the bright fast components (B3 and B6, respectively). 
Figure 44 shows knot D ~6 mas from the core. This can be associated with a feature 
(A PT ) imaged at longer wavelengths by Pauliny-Toth (1998). According to Pauliny-Toth 
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(1998) Apt accelerated in 1984-1989 from 16c to 30c and then decelerated to 20c in 1992. 
Our observations suggest a possible continuation of the deceleration such that the apparent 
speed of the knot has decreased to 6c. 
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Table 1. Sample and Parameters of Maps shown in Figures 1-15 
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Table 2. 


Log of VLBA Observations 
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Table 3. Best-Fit Polynomials for Component CI of 3C 66A 
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Table 4. Best-fit Polynomials for Component BI of 3C 273 
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Table 5. Jet Velocities 
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.84±0.06 


B5 


8 


0. 


.16 





.94 


57.2±1.6 


1 


56, 


1±3. 


.1 


0, 


98±0. 


.06 


3.3±0.2 


1999. 


,23±0.03 


54 


7 


0. 


.25 


0. 


66 


59.0±0.8 


1 


58, 


,8±5, 


.2 


1. 


,00±0, 


.10 


3.4±0.3 


1999 


,05±0.03 


B'S 


8 


0. 


.20 





.93 


63.0±1.1 


1 


62, 


4±3. 


.9 


1. 


.39±0. 


.08 


4.7±0.3 


1998. 


.76±0.03 


B2 


6 


0. 


.03 


0. 


.91 


64.4±0.9 


1 


65, 


5±6. 


.0 


1. 


,11±0. 


.14 


3.7±0.5 


1997 


.85±0.05 


Bl 


4 





.05 


1. 


.08 


58.9±0.9 


1 


62.6±12, 


.3 


0. 


83±0. 


.26 


2.8±0.9 


1997 


.24±0.11 


CI 


13 


0. 


.20 


4. 


.38 


65.8±1.5 


1 


60 


9±2 


.9 


1. 


.77±0. 


.06 


6.0±0.2 


1996. 


,93±0.24 


cl 


16 


0. 


.09 


3. 


.73 


67.1±2.7 


2 


56, 


,9±0. 


.9 


1. 


07±0, 


.02 


3.6±0.1 


1996 


.31±0.62 


c2 


9 


0. 


.07 


5. 


.29 


65.3±1.8 


1 


54, 


,5±2. 


.3 


1 


23±0. 


.04 


4.1±0.2 


1996 


,1±0.9 


Al 


17 


0. 


.58 


0. 


.12 


63.5±4.7 





76. 


,8±0. 


.7 


0. 


02±0. 


.03 


0.05±0.09 






A2 


12 





.37 


0. 


.27 


60.2±3.5 





31 


,8±1. 


.8 


0. 


02±0, 


03 


0.07±0.10 






B2 


4 


1. 


.30 





.13 


-79.2±8.1 


1 


-71. 


,7±5. 


.2 





16±0. 


.07 


7.6±3.5 


2000 


,22±0.35 


Bl 


14 


0. 


.27 





32 


-147±25 


3 


171, 


,8±2. 


.2 


0. 


22±0. 


.06 


10.9±2.8 


1998. 


.59±0.16 


V 


4 


0. 


.09 


1 


.03 


-117.8±1.4 


1 


-116. 


,1±8. 


.9 


1 


96±0. 


.20 


4.3±0.4 


2000 


,36±0.06 


u 


8 


0. 


.27 


1 


.28 


-123.1±1.9 


1 


-122. 


,9±0. 


.5 


1 


69±0. 


.07 


3.7±0.2 


1999. 


,70±0.02 


ul 


5 





.13 


1. 


.10 


-117.7±4.6 


1* 


-120, 


,2±7, 


.5 


1. 


.84±0. 


.15 


4.0±0.3 


1999 


,92±0.02 


t 


10 


0. 


.17 


1 


.95 


-125.1±5.9 


1 


-128. 


,6±5. 


.2 


2 


04±0. 


.08 


4.5±0.2 


1999 


,27±0.03 


ol 


16 


0. 


.27 


3. 


.28 


-121.4±3.0 


2 


-118, 


,6±2. 


.8 


2. 


03±0, 


.06 


4.5±0.1 


1998 


.07±0.03 


o2 


14 


0. 


.42 


1. 


.85 


-122.9±3.2 


2 


-119. 


,1±1, 


.4 


1 


,52±0. 


03 


3.4±0.1 


1998. 


,14±0.07 


I 


7 


0. 


.16 


2 


.85 


-124.9±1.9 


1 


-120. 


,9±6. 


.8 


1 


66±0. 


13 


3.6±0.3 


1997. 


,03±0.14 


h 


16 


0. 


.16 


5. 


.49 


-113.2±4.1 


1* 


-103 


,8±1, 


.7 


1. 


82±0. 


.06 


4.0±0.1 


1996, 


.35±0.19 


d 


7 


0. 


.06 


6 


.95 


-113.0±2.0 


1 


-95. 


,0±5. 


.9 


1 


66±0. 


.45 


3.6±1.0 


1994. 


.88±0.28 


Al 


15 





.22 





.19 


-122.6±10.9 


1 


34. 


3±1. 


.7 


0. 


07±0. 


.03 


0.15±0.06 






B8 


11 


1 


33 





.13 


77±16 


2 


44. 


3±1. 


.0 


0. 


11±0. 


.02 


9.3±1.4 


1999 


,12±0.50 


B8p 


9 


0. 


.04 





.19 


60±22 


1 


33 


,6±2. 


.5 


0. 


22±0. 


.04 


18.6±3.7 


1999 


.61±0.23 
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Table 5 — Continued 



Source 


Knot N <S> < R > 


< e >° 


i 


< $° > 


< ll > < (3 app > 




Jy mas 








mas yr _1 c 



OJ 287 



3C 273 



3C 279 



1510-089 



Bl 


6 


0. 


.43 


0. 


.16 


108±22 


1* 


66.9±3. 


5 





,25±0 


.08 


20 


,8±7 


,1 


1998.69±0.55 


B6 


13 


0. 


.32 


0. 


,34 


86±20 


3 


52.9±3. 


2 





,26±0, 


.04 


21 


,6±3 


.1 


1997.95±0.28 


B5 


10 


0, 


.08 


0. 


,43 


75±12 


1* 


53.1±2, 


.0 





,25±0 


.04 


20 


,7±3 


.1 


1997.44±0.62 


BA 


9 


0. 


.71 


0, 


.39 


69±14 


3 


36.6±5, 


,6 





,21±0 


.08 


17. 


,8±8. 


,2 


1995.58±0.15 


B3 


6 


0. 


.65 





,53 


68±17 


1* 


31.9±2, 


.5 





,31±0 


.04 


26 


,1±3 


.7 


1995.46±0.85 


B6 


4 


0. 


,28 





.39 


-113.9±11.0 


1* 


-131.9±12. 


.5 





,93±0 


.22 


18. 


,0±4 


.1 


1999.57±0.11 


B5 


6 


0. 


,28 


0, 


,39 


-103.6±22.1 


1* 


-117. 6±4, 


.5 





,60±0 


.10 


11 


,6±1 


.9 


1998.61±0.09 


BA 


5 


0. 


,27 


0, 


,33 


-99.2±20.5 


1* 


-115. 9±6. 


.5 





,92±0 


.14 


17 


,8±2 


.8 


1998.23±0.04 


Al 


17 


0. 


,24 


0. 


,10 


-90±27 


1* 


-153.4±0. 


,4 





,03±0 


,01 


0, 


,5±0, 


,2 


1992.67±5.0 


A2 


17 


0. 


,41 


0, 


,30 


-111±4 


1* 


-101.4±0. 


,4 





,03±0 


.02 


0. 


,6±0. 


.3 


1989.7±5.0 


A2* 


24 


0, 


,36 


0, 


,28 


-109±9 


1* 


-178. 2±6, 


.5 





,02±0 


.01 


0, 


,4±0. 


.2 


-1970.5 


A3 


16 


0. 


,12 


1. 


,00 


-108.9±2.2 


1 


69.3±1, 


,2 





,06±0 


.03 


1. 


,3±0, 


.6 




B2 


10 


3 


,29 


1 


,16 


-129.6±4.5 


2 


-123. 5±3, 


3 


1 


,18±0 


,06 


12 


,2±0 


.6 


1999.25±0.02 


62 


6 


1. 


,70 


1, 


,19 


-121.1±1.0 


1* 


-127.8±8. 


.1 





,36±0 


.15 


3 


,7±1. 


.5 


1997.5±0.4 


Bl 


16 


7, 


,58 


1, 


.61 


-127.7±4.0 


4 


-120.3±3. 


.0 





,79±0 


.06 


8, 


,2±0, 


.6 


1997.89±0.16 


bl 


17 


4, 


,04 


1 


,16 


-122.6±5.0 


1 


-119.4±1. 


.5 





,47±0 


.03 


4, 


,9±0, 


.3 


1997.27±0.16 


b s 


8 


0. 


,55 


2. 


,67 


-135.1±2.4 


1 


-139. 6±2. 


.8 


1 


,12±0 


.05 


11 


,6±0 


.5 


1998.13±0.18 


Al 


17 


1. 


,14 


0. 


,16 


-131.8±12.3 





-11.8±0. 


.6 





,00±0 


.02 


0, 


,0±0. 


.2 




A2 


11 


5. 


,74 





,78 


-123.3±3.8 


1* 


-104.1±1. 


.7 





,06±0 


.06 


0. 


,6±0, 


,6 


1984±12 


CI 


10 


0. 


,32 


8. 


,27 


-114.5±1.2 


1* 


-123. 9±8. 


.9 





,67±0 


.16 


6, 


,9±1, 


.7 


1988. 2±2. 2 


C16 




3 


33 


0. 


,30 


— 144.9±7.3 


1 


_1 QC Q-LR 

lOU. O —1— \J ■ 


.8 


o 


,53±0 


.11 


16 


.9±3 


.5 


2000 27+0 05 


C15 


8 


2. 


,51 


0. 


,42 


-137.4±4.9 


1 


-132.4±4. 


.3 





,54±0 


.07 


17. 


,2±2 


.3 


1999.85±0.04 


C14 


10 


1. 


,34 


0. 


,53 


-134.5±5.5 


1 


-128. 8±4, 


.2 





,54±0 


.08 


17 


,2±2 


.4 


1999.50±0.09 


C13 


13 


2. 


,03 


0, 


,59 


-129.8±4.2 


1 


-126. 5±2, 


,9 





,51±0 


,05 


16 


,2±1. 


,6 


1998.99±0.07 


C12 


11 


2, 


,02 


0. 


,54 


-129.5±3.1 


1 


-124.2±2. 


3 





,49±0 


.04 


15 


,7±1 


.4 


1998.50±0.08 


Cll 


7 


1. 


93 


0, 


,38 


-135.0±3.6 


1 


-129.4±5. 


.8 





,32±0 


.10 


10 


,1±3. 


.2 


1997.59±0.12 


CIO 


14 


1. 


.75 





,70 


-131.1±5.9 


1 


-119. 2±1. 


.6 





,35±0, 


.03 


11 


,1±1. 


.1 


1997.23±0.47 


C9 


21 


2. 


.23 


0, 


,98 


-129.8±4.9 


1 


-123.4±0. 


,4 





,40±0 


.01 


12 


,8±0 


.3 


1996.89±0.22 


C8* 


27 


3 


.13 


0, 


,81 


-129.3±3.2 


3 


-126. 6±0. 


,7 





,33±0 


.01 


10 


,5±0 


.4 


1996.10±0.31 


CA 


46 


1. 


.28 


3 


,54 


-116.2±2.0 


3 


-127.9±0. 


.8 





,31±0 


.01 


9. 


,9±0. 


.5 


1987.7±2.7 


CAa 


5 


1 


.33 


2, 


.66 


-115.0±1.0.3 


1 


-97.5±2. 


1 





,41±0 


.13 


13 


,0±4 


,2 


1988.2±9.8 


CAb 


11 


1. 


.56 


3, 


,19 


-114.3±0.5 


1 


-123. 8±3. 


.7 





18±0. 


.07 


5, 


,7±2. 


,2 


1979.1±3.4 


CAc 


30 


1. 


.26 


3 


.81 


-117.0±1.8 


1 


-140. 2±1. 


.3 





,35±0 


,02 


11 


,1±0 


.7 


1987.5±4.5 


E 


4 


1. 


.12 


0. 


,28 


-136.2±4.2 


1* 


-147.8±1. 


.0 





,04±0, 


.02 


1, 


,3±0. 


,5 


1992.88±1.41 


53 


5 


0. 


.24 


0. 


,55 


-46.6±6.4 


1 


-41.9±6. 


.0 


1 


,25±0 


.10 


28 


,0±2 


.4 


1999.67±0.03 


B2 


4 


0. 


.10 


0. 


,64 


-33.8±2.3 


1 


-29.5±7. 





1 


,56±0 


.13 


35 


,0±3 


,0 


1998.81±0.04 
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Table 5 — Continued 



Source 


Knot N < 5 > < R > 


< e >° 


i < $° > 


< n > < (3 app > 


To 




Jy mas 






mas yr _1 c 





3C 345 



1803+784 



1823+568 



BL Lac 



CTA 102 



3C 454.3 



Bl 


4 


0. 


.04 


0. 


.78 


-48.5+3. 


.6 


1 


-43 


.2+9.3 


2. 


.04+0, 


,16 


45. 


,9+3 


.6 


1998. 


.11+0.09 


Al 


15 


0. 


.28 


0. 


.14 


-49.7+21. 


,6 


1* 


-72, 


.7+2.6 


0. 


,02+0, 


,02 


0. 


,4+0 


.4 






A2 


16 





.05 


0. 


.54 


-33.9+7. 


.6 


1* 


105, 


.8+0.9 


0, 


,03+0. 


.02 


0. 


,5+0 


.4 






5 


9 


0. 


.18 





,13 


127.1+32. 


.8 


1* 




17+1 


0. 


.07+0. 


.02 


1, 


.6+0 


.5 






CIS 


5 





.96 


0. 


.40 


-87.7+2. 


.3 


1* 


-84, 


.9+1.3 


0, 


.42+0. 


.09 


14. 


,6+3 


.3 


2000. 


,00+0.34 


C12 


7 


0. 


.74 





,59 


-86.5+2. 


.3 


1* 


-88. 


.6+0.2 





.66+0. 


.05 


23 


1+1 


.7 


1999. 


,80+0.59 


Cll 


10 


0. 


.78 


0. 


.69 


-95.2+2. 


.7 


2 


-94, 


.5+0.4 


0. 


.69+0, 


.03 


23 


,9+1 


.1 


1999. 


.24+0.29 


CIO 


14 


0. 


.35 


0. 


.83 


-95.0+4. 


.6 


1 


-102, 


.8+0.6 


0, 


.55+0. 


.02 


19. 


,2+0 


.7 


1998 


.60+0.41 


C9 


17 


0. 


.40 


1 


.23 


-106.2+7. 


.5 


4 


-79. 


.1+1.0 





,45+0. 


.04 


15. 


,6+1 


.4 


1996 


,71+0.26 


c9 


8 


0. 


.78 


1 


.23 


-110.4+2. 


.9 


1 


-119. 


.6+2.1 


0. 


.27+0. 


.04 


9. 


.2+1 


.5 


1996 


.11+0.57 


C8 


10 


0. 


.33 


1. 


.54 


-97.2+5. 


,6 


3 


-65. 


.6+4.7 


0. 


.32+0. 


.10 


11 


.0+3 


.6 


1994. 


.0+0.8 


C7 


16 


0. 


.14 


2. 


.89 


-96.8+1. 


,4 


2 


-81, 


.4+0.3 


0. 


,29+0, 


.02 


10. 


.1+0 


.5 


1993 


.77+0.36 


Al 


14 


2. 


.29 


0. 


.10 


-84.6+5. 


,2 





-100, 


.4+0.5 


0. 


,02+0, 


.02 


0. 


.7+0 


.6 






B2 


13 


0. 


.12 





.38 


-62.2+20. 


,8 


4 


-75. 


.4+1.6 


0. 


.42+0. 


.05 


15. 


.9+1 


.9 


1998. 


.58+0.15 


Bl 


14 


0. 


.14 


0. 


.70 


-77.5+5. 


,1 


4 


-60. 


.0+2.4 


0, 


.42+0. 


.05 


15. 


,7+1 


.8 


1997 


.64+0.08 


Al 


17 


0. 


.21 


0. 


,28 


-77.9+4. 


.5 


1* 


150. 


.0+0.7 


0. 


,01+0. 


.01 


0, 


,4+0 


.5 






A2 


17 





.11 


1 


,44 


-90.1+2. 


.2 


1* 


-36 


.1+0.8 


0, 


.04+0. 


.01 


1 


,5+0 


.5 






CI 


7 


0. 


.12 


1. 


,66 


-96.0+1. 


,1 


1 


-93, 


.2+0.3 


0, 


.28+0. 


.04 


10, 


,5+1 


.5 


1992. 


.6+3.6 


Al 


17 


0. 


.54 


0. 


,06 


-163.2+10 


.4 





-23. 


.4+0.4 





,00+0. 


.01 


0, 


.1+0 


.3 






A2 


17 


0. 


.07 


0. 


,21 


-164.4+5. 


,2 





-40, 


.9+0.5 


0. 


.02+0. 


.01 


0, 


,6+0 


.3 






A3 


17 


0. 


.11 





,52 


1 R2 5+2 


2 


o 


—89. 


.7+0.6 


0, 


.02+0. 


,01 


0. 


.7+0 








C2 


11 


0. 


.04 


1 


.61 


-162.1+1. 


,1 


1 


-174. 


.9+0.3 





.20+0. 


.02 


7, 


.7+0 


.6 


1991 


.5+1.0 


CI 


16 





.04 


1. 


.91 


-160.8+2. 


.4 


1 


-153, 


.1+0.6 


0. 


.10+0. 


.01 


3. 


.8+0 


.5 


1979. 


.4+1.3 


513 


10 





.10 


2 


.27 


-158.3+5. 


,1 


3 


-170. 


.2+2.9 


1 


.95+0. 


.11 


9, 


,0+0 


.5 


1999 


.29+0.04 


512 


9 


0. 


.10 


1. 


.83 


-158.3+1. 


,2 


1 


-163, 


.2+2.4 


1. 


.82+0. 


.06 


8, 


,4+0 


.3 


1998, 


.93+0.12 


511 


10 


0. 


.11 


1 


.58 


-159.3+4. 


.4 


3 


-167. 


.1+4.2 


1 


.36+0. 


.13 


6. 


,3+0 


.6 


1998. 


.02+0.12 


510 


10 


0. 


33 


2 


.20 


-167.6+6. 


,1 


2 


-175. 


.8+2.4 


0. 


.91+0. 


03 


4, 


,2+0 


.2 


1996 


,93+0.08 


Al 


17 





.73 


0. 


.10 


-160.4+9. 


.7 





-178, 


.9+0.3 


0. 


.01+0 


.01 


0, 


.1+0 


.1 






A2 


17 


0. 


.34 


0. 


.29 


-158.9+6. 


.4 


1* 


143. 


.9+0.6 





.01+0. 


.01 


0. 


,1+0 


.1 






B6 


10 


0. 


.92 


0. 


.30 


118.2+9. 


.4 


1 


121, 


.9+1.9 


0. 


.34+0. 


.04 


18. 


.2+1 


.9 


1999. 


.54+0.04 


b6 


5 


0. 


.24 





.32 


121.1+2. 


.4 


1 


125. 


.0+6.0 


0. 


.16+0. 


,11 


8. 


.5+6 


.0 


1998. 


.9+1.5 


B5 


16 


0. 


.99 


0. 


.50 


112.8+3. 


.9 


1 


119. 


.7+0.9 





.29+0. 


,02 


15. 


.4+0 


.9 


1997. 


.9+0.2 


bb 


12 


0. 


.29 


0. 


.52 


115.3+3. 


.9 


1* 


123, 


.8+1.4 


0. 


.23+0. 


.03 


12. 


.5+1 


.5 


1998, 


.0+0.3 


Al 


16 


0. 


.59 


0. 


.12 


117.1+13. 


.0 





-129. 


.5+0.7 


0. 


.01+0. 


.01 


0. 


.6+0 


.7 






C 


15 


0. 


.13 


1. 


.99 


137.6+1. 


.6 





96, 


.2+0.7 


0. 


.01+0 


.05 


0. 


.6+2 


.6 






B6 


7 


1. 


.64 





.31 


-101.0+21 


,8 


2 


-79. 


.5+1.2 





.53+0. 


.05 


24. 


.8+2 


.5 


1999 


.80+0.37 
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Table 5 — Continued 



Source Knot 


N 


< S > 


<R> 


< e >° 


I 


< $° > 


< n > 


^ ftapp ^ 


T 






Jy 


mas 








mas yr -1 


C 




66 


4 


0.65 


0.26 


-92.1±2.8 





-67.2±3.8 


0.03±0.09 


1.5±4.3 




Bb 


5 


0.31 


0.19 


-85.1±4.7 


1 


-88.7±0.3 


0.22±0.09 


10.0±4.3 


1999.04±0.44 


BA 


15 


0.68 


0.61 


-86.2±10.5 


2 


-98.3±0.4 


0.41±0.02 


19. Oil. 1 


1998.36±0.07 


B3 


7 


0.85 


1.37 


-74.6±1.2 


1 


-78.1±2.7 


0.28±0.11 


13.2±5.3 


1995.70±0.27 


63 


6 


0.79 


0.73 


-90.9±2.7 


1* 


-102.1±3.5 


0.23±0.15 


10.7±6.8 


1995. 3±2. 2 


C 


17 


1.27 


0.63 


-79.5±4.1 


1* 


2.7±0.2 


0.01±0.01 


0.4±0.6 




c* 


29 


1.44 


0.63 


-73.4±8.4 


1* 


-165.1±0.3 


0.03±0.01 


1.6±0.3 




N 


14 


0.37 


0.86 


-43.9±7.0 


1* 


-14.6±0.9 


0.16±0.02 


7.7±1.0 


1992±5 


D 


16 


0.11 


5.80 


-83.1±1.3 


1 


-21.6±1.9 


0.12±0.04 


5.5±1.9 
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Table 6. Acceleration in the Jets 



Source Knot N < \i > < <J>° > fi\\ fi ± 

mas yr -1 mas yr -2 mas yr 



0420-014 


Bl 


14 





.22+0. 


.06 


171. 


.8+2. 


.2 


0. 


.07+0, 


.02 


0. 


.02+0.01 


3C 111 


cl 


16 


1 


.07+0. 


.02 


56. 


,9+0. 


.9 


-0 


.22+0, 


.02 


-0. 


.26+0.02 


3C 120 


ol 


16 


2 


.04+0. 


.06 


-118 


.6+2 


.8 


0. 


.29+0, 


.02 


0. 


,16+0.02 




o2 


14 


1 


.52+0. 


.03 


-119 


.1+1 


.4 


-0 


,66+0, 


.03 


-0. 


.12+0.04 


0528+134 


B8 


11 





.11+0. 


.02 


44. 


.3+1. 


.0 


0. 


,10+0, 


.01 


-0. 


.07+0.02 




B6 


13 





.26+0. 


.04 


52. 


.9+2. 


.2 


0, 


,06+0, 


.03 


-0. 


.17+0.04 




B4 


9 





.21+0. 


.08 


36 


,6+5. 


.6 





,11+0. 


.02 


-0 


.09+0.03 


3C 273 


Bl 


16 





.79+0. 


.06 


-120. 


.3+3 


.0 


-0, 


.14+0. 


.03 


0. 


,10+0.02 




B2 


10 


1 


.18+0. 


.06 


-123 


.5+3 


.3 


0. 


.04+0. 


.03 


0. 


.26+0.04 


3C 279 


C8* 


27 





.33+0. 


.01 


-126 


.6+0, 


.7 


0, 


.07+0. 


.01 





.02+0.01 




C4 


46 





.31+0. 


.01 


-127. 


.9+0 


.8 


0. 


.02+0. 


.01 


-0 


.05+0.01 


3C 345 


Cll 


10 





.69+0. 


.03 


-94. 


,5+0. 


.4 


0, 


.22+0. 


.03 


0. 


.19+0.04 




C9 


17 





.45+0. 


.04 


-79. 


.1+1. 


.0 


0. 


.11+0. 


.02 


-0 


.08+0.03 




C8 


10 





.32+0. 


.10 


-65. 


.6+4. 


.7 


-0. 


.20+0. 


.03 


-0 


.23+0.03 




Cl 


16 





.29+0. 


.02 


-81. 


.4+0. 


.3 


-0. 


.07+0. 


.01 


-0. 


.01+0.01 


1803+784 


B2 


13 





.42+0. 


.05 


-75. 


.4+1. 


6 


0. 


.30+0. 


.05 


-0 


.10+0.04 




Bl 


14 





.42+0. 


.05 


-60. 


.0+2. 


.4 


0. 


,52+0. 


.04 


-0. 


.29+0.03 


BL Lac 


513 


10 


1 


.95+0. 


.11 


-170 


.2+2 


.9 


-0. 


.21+0, 


.07 


-0. 


.76+0.06 




511 


10 


1 


.36+0. 


.13 


-167 


.1+4. 


.2 


0. 


,07+0, 


.06 


-0 


.62+0.05 




510 


10 





.91+0. 


.03 


-175 


.8+2, 


.4 


-0. 


.32+0, 


.03 


-0. 


.71+0.04 


3C 454.3 


B6 


7 





.53+0. 


.05 


-79. 


,5+1. 


.2 


0. 


,37+0, 


.06 


0, 


,24+0.07 




B4 


15 





.41+0. 


.02 


-98. 


.3+0, 


.4 


-0 


.04+0, 


.01 


-0. 


.09+0.02 
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Table 7. Parameters of Total and Polarized Intensity Maps shown in Figure 17 



Source 


Beam 


P.A. 


-^pcak 


T p 

peak 




rms(I) 


rms(J p ' 




mas x mas 





mjy Beam -1 


mjy Beam -1 


mjy Beam -1 


mjy Beam -1 


mjy Bean 


3C 66A 


0.32x0.19 


-13 


500 


30 


1.5 


0.7 


0.5 


3C 111 


0.30x0.18 


-10 


803 


6 


2.0 


0.9 


0.9 


0420-014 


0.36x0.14 


-8 


3504 


112 


10.5 


4.0 


2.7 


3C 120 


0.38x0.15 


-6 


650 


11 


2.9 


0.7 


1.4 


0528+134 


0.30x0.15 


-7 


1426 


53 


3.8 


1.8 


1.2 


OJ 287 


0.30x0.14 


-8 


1115 


32 


4.4 


2.5 


2.6 


3C 273 


0.42x0.17 


-6 


7178 


655 


35.9 


7.6 


14.1 


3C 279 


0.38x0.14 


-7 


9695 


432 


24.2 


3.5 


3.1 


1510-089 


0.40x0.15 


-5 


1517 


36 


7.6 


1.5 


3.4 


3C 345 


0.22x0.15 


-10 


1470 


144 


7.4 


1.2 


1.4 


1803+784 


0.16x0.15 


30 


1211 


61 


4.2 


1.3 


2.2 


1823+568 


0.19x0.15 


-11 


781 


55 


3.9 


0.6 


1.4 


BL Lac 


0.28x0.15 


-4 


2709 


128 


6.8 


2.2 


2.7 


CTA 102 


0.24x0.16 


-9 


2456 


109 


6.1 


3.1 


4.6 


3C 454.3 


0.30x0.15 


-6 


3729 


206 


9.3 


1.8 


2.9 
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Table 8. Jet Components in Quasars 



Source Knot r var 9™(°) «5 var T b , obs (10 10 if) T hMt (10 9 K) 



0420-014 


B2 


10.4 


2.4 


<17.5 


4.69 


1.09 




Bl 


11.5 


3.6 


<15.0 


1.64 


0.49 


0528+134 


B8 


27.4 


0.4 


<53.2 


7.94 


0.62 




Bl 


35.2 


0.5 


<63.6 


2.17 


0.125 




B6 


25.5 


1.2 


39.2 


2.58 


0.34 




B5 


21.1 


3.3 


17.3 


0.38 


0.20 




BA 


27.5 


0.8 


48.5 


3.72 


0.34 




B2> 


39.5 


0.6 


69.2 


1.07 


0.05 


3C 273 


B2 


13.8 


6.9 


7.4 


15.4 


6.62 




Bl 


8.3 


6.1 


9.4 


16.4 


4.66 




CI 


8.2 


3.8 


12.6 


0.23 


0.04 


3C 279 


C16 


17.1 


2.9 


19.7 


20.2 


2.65 




C15 


18.2 


2.3 


23.9 


25.3 


2.38 




CIA 


17.4 


2.8 


20.0 


9.37 


1.19 




C13 


16.3 


3.3 


17.6 


16.1 


2.57 




C12 


16.6 


2.5 


22.0 


23.3 


2.52 




Cll 


12.9 


2.2 


<20.9 


12.6 


1.49 




CIO 


14.9 


1.7 


24.7 


4.96 


0.44 




C9 


16.7 


1.6 


27.3 


3.74 


0.28 




C8 


17.2 


0.6 


33.4 


6.47 


0.35 




C4 


10.5 


3.9 


13.8 


1.05 


0.25 


1510—089 


B3 


30.1 


1.3 


41.3 


1.10 


0.03 






IK 
oO.Z 


1.0 


oy.U 


n 73 
u. / o 






Bl 


47.9 


1.6 


34.2 


0.21 


0.01 


3C 345 


C13 


16.2 


2.3 


23.1 


7.41 


0.79 




C12 


24.6 


1.6 


32.9 


3.25 


0.19 




Cll 


29.8 


3.8 


12.0 


7.37 


2.39 




C10 


19.3 


3.3 


17.6 


12.1 


2.05 




C9 


16.2 


4.6 


11.9 


1.51 


0.49 




C8 


19.9 


0.9 


36.5 


0.50 


0.025 




C7 


13.3 


2.0 


21.9 


0.02 


0.003 


CTA 102 


B6 


20.5 


1.7 


<29.8 


3.54 


0.37 




B5 


15.7 


3.0 


18.8 


14.4 


3.29 


3C 454.3 


B6 


25.2 


1.9 


29.7 


17.6 


1.58 




55 


14.8 


1.5 


<25.8 


1.58 


0.18 




BA 


20.0 


3.9 


13.9 


4.54 


1.49 




B'S 


13.4 


3.7 


15.4 


24.7 


6.76 
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Table 8 — Continued 



Source Knot 


pvar 


err) 


<5 var 


T b , obs (10 10 K) 


T b ant(10 9 K) 


C 


14.2 


0.2 


28.2 


6.76 


0.66 


D 


10.2 


1.7 


18.7 


0.03 


0.007 



Table 9. Jet Components in BL Lac Objects 



Source 


Knot 


pvar 


er(°) 


^var 


T b ,obs(10 10 K) 


T b ,i„t(10 9 K) 


3C 66A 


CA 


5.0 


1.8 


<9.7 


0.73 


0.285 




C3 


5.9 


2.8 


<10.9 


0.44 


0.141 




C2 


7.7 


2.8 


<13.5 


0.35 


0.079 




CI 


11.4 


0.8 


22.2 


0.15 


0.014 




B7 


29.7 


3.3 


<15.4 


0.21 


0.037 




B6 


26.1 


4.0 


12.1 


0.12 


0.032 




B4 


27.6 


2.5 


22.1 


0.02 


0.002 


OJ 287 


B6 


20.7 


4.7 


10.6 


0.41 


0.117 




B5 


12.3 


3.4 


16.1 


1.31 


0.183 




54 


19.9 


1.8 


<28.8 


1.70 


0.088 


1803+784 


B2 


17.2 


2.2 


23.7 


0.34 


0.038 




Bl 


15.7 


3.6 


16.1 


0.25 


0.054 




CI 


20.8 


0.8 


38.7 


0.02 


0.001 


1823+568 


C2 


9.8 


2.9 


15.7 


0.20 


0.005 




CI 


13.8 


0.6 


27.0 


0.14 


0.012 


BL Lac 


513 


9.1 


7.0 


8.1 


0.03 


0.009 




512 


8.6 


7.9 


7.2 


0.16 


0.061 




511 


7.2 


13.2 


3.9 


0.07 


0.082 




510 


5.1 


6.0 


8.0 


0.07 


0.025 



Table 10. Jet Components in Radio Galaxies 



Source Knot r var Q™(°) T b , obs (10 10 K) T b a nt (10 9 K) 



B7 


2.8 


24.8 


2.3 


1.49 


3.82 


B6 


3.7 


22.9 


2.4 


8.00 


0.20 


Bb 


3.5 


19.2 


<3.0 


1.44 


2.40 


BA 


4.0 


24.4 


2.1 


1.61 


4.94 


S3 


4.8 


10.4 


<5.5 


0.06 


0.04 


B2 


4.2 


21.6 


2.5 


0.29 


0.69 


Bl 


3.2 


12.5 


<4.3 


0.23 


0.21 


CI 


6.6 


13.0 


4.1 


0.15 


0.15 


V 


6.6 


23.0 


1.7 


0.19 


0.85 


u 


4.3 


22.3 


2.3 


0.75 


1.89 


t 


4.9 


16.9 


3.2 


0.11 


0.16 


ol 


7.0 


22.0 


1.7 


3.69 


15.2 


I 


4.0 


21.4 


2.5 


0.88 


1.91 


h 


4.5 


19.8 


2.7 


0.31 


0.61 


d 


3.9 


11.3 


4.9 


0.03 


0.02 
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Table 11. Global Parameters of the Jets 



Source 


D(Gpc) 








< e > 


O 




6° 




< r > 




<5> 




/ej 


3C 66A 


2.46 


-162.0+25 


9.5+4.6 


2, 


,1+1. 


.0 





.4+0, 


,2 


10 


.4+3. 


5 


16.0+4. 


.9 


1.9 


3C 111 


0.222 


63.5+4. 


5 


9.1+2.7 


18 


.1+5. 


.0 


2 


.8+0 


.8 


4, 


.4+1. 


3 


3.4+1. 


.1 


2.4 


0420-014 


5.925 


-139.2+21. 


.9 


21+10 


3. 


.0+0. 


.6 


1 


.1+0 


.5 


11 


.0+0. 


5 


16.2+1. 


.2 


1.1 


3C 120 


0.145 


-121.7+6. 


.1 


10.9+3.8 


20 


.5+1 


.8 


3 


.8+1, 


.3 


5. 


.3+1. 


2 


2.4+0. 


.6 


2.4 


0528+134 


16.12 


74.7+22. 


.4 


35.1+9.3 


1, 


.1+1. 


.0 





.7+0 


,2 


28 


.3+4. 


5 


47.5+13. 


.8 


7.0 


OJ 287 


1.59 


-105.8+16. 


.9 


11.1+5.4 


3 


.2+0. 


.9 





.8+0 


.4 


16 


.5+4. 





18.9+6. 


.4 


1.3 


3C 273 


0.76 


-126.4+8. 


.6 


13.4+2.5 


6. 


,1+0. 


.8 


1 


.4+0. 


.3 


10 


.6+2. 


.8 


9.0+1. 


A 


0.7 


3C 279 


3.10 


-131.5+10. 


.8 


8.3+2.7 


2. 


,1+1. 


.1 





.4+0 


,2 


15 


.5+2. 


5 


24.1+6. 


.5 


3.5 


1510-089 


1.93 


-42.6+16. 


.4 


4.8+6.5 


1. 


,4+0. 


.4 





.2+0. 


.2 


36 


.6+7. 





38.6+2. 


,8 


2.6 


3C 345 


3.50 


-94.9+8. 


6 


10.7+4.4 


2. 


,7+0. 


.9 





.5+0 


.2 


18 


.7+5. 


5 


20.2+5. 


.0 


3.7 


1803+784 


4.11 


-78.6+31. 


.1 


7.5+3.6 


2. 


.1+1. 


.2 





.3+0 


.2 


18 


.1+2. 


3 


27.0+9. 


.9 


1.0 


1823+568 


4.00 


-162.6+5. 


.9 


2.4+2.0 


1, 


,6+1. 


.1 





.1+0. 


.2 


12 


.0+2. 





21.9+5. 


6 


0.5 


BL Lac 


0.31 


-161.1+7. 


.3 


13.9+4.3 


7, 


,7+1. 


.9 


1 


.9+0 


,6 


7. 


.0+1. 


.8 


7.2+1. 


.1 


1.1 


CTA 102 


6.92 


123.1+19. 


.9 


4.8+1.9 


2. 


.6+0, 


.5 





.2+0 


.1 


17, 


.2+2. 





22.3+4. 


,5 


1.2 


3C 454.3 


5.48 


-76.8+27. 


.8 


36.8+9.8 


1. 


,3+1, 


,2 





.8+0. 


.2 


15. 


.6+2. 


2 


24.6+4. 


,5 


1.9 
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Fig. 1. — Total intensity images of 3C 66A at 43 GHz. North is toward the top and east 
is toward the left. Line segments within each image indicate direction of electic vectors, 
with length proportional to the polarized intensity. Line segments outside the images give 
direction of polarization at 3 mm (dark gray), 0.85/1.3 mm (bold black), and at optical 
wavelengths (light gray) with length reflecting the percent polarization. Percent polarization 
is printed near corresponding line segments along with the epoch of observation. Darkness 
of the print indicates the wavelength of observation. If epochs of the JCMT and optical 
observations coincide, only the JCMT epoch is indicated. Parameters of the maps are listed 
in Table I. The resolution beam is shown by the cross-hatched ellipse in the lower left 
corner. Designation of knots follows Table 5. The figure can be found at the web-site 
www.bu.edu/blazars/multi.html. 
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Fig. 2. — Left panel: Total intensity images of 3C 111 at 43 GHz along with polarization 
measurements at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelenths. The epochs of 
BIMA observations are marked by the letter B. See caption for Figure 1 for details. Right 
panel: Total intensity images of the inner jet of 3C 111 at 43 GHz. The figure can be found 
at the web-site www.bu.edu/blazars/multi.html. 
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Fig. 3. — Total intensity images of 0420—014 at 43 GHz along with polarization measure- 
ments at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelenths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Fig- 
ure 1. 
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Fig. 4. — Total intensity images images of 3C 120 at 43 GHz along with polarization mea- 
surements at 7 mm, 3 mm, and 1.3/0.85 mm. Measurements of the optical polarization are 
consistent with Galactic interstellar polarization for this sight line and are not included. See 
caption for Figure 1. 
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Fig. 5. — Total intensity images of 0528+134 at 43 GHz along with polarization measure- 
ments at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelenths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Fig- 
ure 1. 
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Fig. 6. — Total intensity images of OJ 287 at 43 GHz along with polarization measurements 
at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelenths. The epochs of BIMA observations 
coincide with the epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 7. — Total intensity images of 3C 273 at 43 GHz along with polarization measurements 
at 7 mm, 3 mm, and 1.3/0.85 mm. The epochs of BIMA observations coincide with the 
epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 8. — Total intensity images of 3C 279 at 43 GHz along with polarization measurements 
at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelenths. The epochs of BIMA observations 
coincide with the epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 9. — Total intensity images of PKS 1510—089 at 43 GHz along with polarization mea- 
surements at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Figure 
1. 
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Fig. 10. — Total intensity images of 3C 345 at 43 GHz along with polarization measurements 
at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA observations 
coincide with the epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 11. — Total intensity images of 1803+784 at 43 GHz along with polarization measure- 
ments at 7 mm, 3 mm, and 1.3/0.85 mm. The epochs of BIMA observations coincide with 
the epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 12. — Total intensity images of 1823+568 at 43 GHz along with polarization mea- 
surements at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Figure 
1. 
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Fig. 13. — Total intensity images of BL Lac at 43 GHz along with polarization measurements 
at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA observations 
coincide with the epochs of corresponding VLB A images. See caption for Figure 1. 
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Fig. 14. — Total intensity images of CTA 102 at 43 GHz along with polarization mea- 
surements at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Figure 
1. 
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Fig. 15. — Total intensity images of 3C 454.3 at 43 GHz along with polarization measure- 
ments at 7 mm, 3 mm, 1.3/0.85 mm, and at optical wavelengths. The epochs of BIMA 
observations coincide with the epochs of corresponding VLBA images. See caption for Fig- 
ure 1. 
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Fig. 16. — Separation of components from the VLBA core vs. epoch of observation. Total 
intensity components are indicated by open circles, polarized intensity components by filled 
circles. In a few cases (3C 66A, 0528+134, 3C 279, and 3C 345) data prior to 1998 are 
from Jorstad et al. (2001a) (triangles), Wehrle et al. (2001) (squares), and Ros et al. (2000) 
(asterisks) (see §A for details). For 3C 279 and 3C 454.3 different symbols are used to 
separate closely located components. The size of a symbol is related to the flux density. 
The solid lines(curves) represent the best polynomial approximation to the data according 
to Table 5. Designation of knots follows Table 5. 
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Fig. 16. — Continued 
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Fig. 16. — Continued 
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Fig. 17. — Total (contours) and polarized (color scale) intensity images at 43 GHz. The 
yellow line segments show the direction of the electric vectors. Parameters of the maps are 
given in Table 7. Designation of components follows Table 5 and Figure 16. The figure can 
be found at the web-site www.bu.edu/blazars/multi.html. 
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Fig. 18. — Trajectories of components that exhibit a statistically significant change of appar- 
ent speed (Table 6). The symbols show measured positions of components while the curves 
represent polynomial approximations to the trajectories. Ballistic trajectories are shown for 
3C 120 (components I and t), OJ 287, and CTA102 (see §A for explanation.) 
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Fig. 19. — Light curves of superluminal jet components. Flux densities are normalized by 
the average flux density listed in Table 5. For each source different components are denoted 
by different symbols. 
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Fig. 19. — Continued 
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Fig. 20. — Size variability timescale At a obs vs. flux variability timescale At var obs for super- 
luminal components in the quasars (open circles), BL Lac objects (filled circles), and radio 
galaxies (triangle). These symbols are applied for all subsequent figures where radio galaxies, 
BL Lacs, and quasars are plotted. The solid line shows the expected relationship during the 
"adiabatic" stage defined by Marscher & Gear (1985). 
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Fig. 21. — Doppler factor 5p estimated from the apparent speed for components with the 
highest apparent speed in the jets vs. Doppler factor 5 var computed from the flux variability 
and size of the components. The solid line shows the relation 6p = 5 var , the dotted line 
shows the least-square fit of the observed dependence. 
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Fig. 22. — Projected transverse size of the jet vs. longitudinal distance from the core. The 
solid lines repesent the best linear fit to the data. 
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Fig. 23. — Apparent speed vs. variability Doppler factor for superluminal components. The 
solid line denotes <5 var = /9 app . The dotted curve encloses the expected location of points 
corresponding to T < 25. 
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Fig. 24. — Distribution of Lorentz factors of superluminal jet components in the 8 quasars, 
5 BL Lac objects, and 2 radio galaxies. 
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Fig. 25. — Inverse Lorentz factor vs. viewing angle for superluminal jet components. 
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Fig. 26. — Left panel: Scatter in the projected jet position angle vs. the average viewing 
angle of jets. Right panel: Projected half opening angle vs. average viewing angle. 
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Fig. 27. — Half opening angle vs. Lorentz factor of the jets. The solid curve corresponds 
to the best approximation according to the x 2 test of the observed dependence within the 
assumed law 9 = p/T, where p is a constant. The dotted curves represent relationships 
between the half opening angle and Lorentz factors expected in the gas-dynamic model for 
different values of the parameter £. 
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Fig. 28. — Left panel: Distribution of intrinsic brightness temperature of the superluminal 
jet components in 8 quasars, 5 BL Lac objects, and 2 radio galaxies. Right panel: Intrinsic 
brightness temperature vs. deprojected distance from the core. 
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Fig. 29. — Left panel: Average Doppler factor vs. average viewing angle of the jets. Right 
panel: Average viewing angle vs. intrinsic half opening angle of the jets. In both panels the 
bold crosses indicate average parameters and their standard deviations for each subclass. 
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Fig. 30. — Apparent speed vs. distance from the core, as defined by a polynomial approxi- 
mation, for accelerating/decelerating knots (Table 6). Short arrows represent the direction 
of the apparent velocity, long arrows indicate the local jet direction; both are derived from 
the best-fit polynomial. Error bars for the apparent speed are computed using uncertainies 
of polynomial coefficients. 
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Fig. 31. — Distribution of change of the instantaneous apparent speed relative to the average 
apparent speed with distance from the core (deprojected). A value of +1 is assigned if the 
instantaneous apparent speed is higher than the average, and —1 if it is lower (see text). 
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Fig. 32. — Distributions of the flux density of the jet components relative to the core - left 
panel: for fast moving knots; right panel: for slow moving knots. 
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Fig. 33. — Apparent speed vs. distance from the core as defined by a polynomial approxi- 
mation, for trailing components observed at >10 epochs; error bars for the apparent speed 
are computed using uncertainies of polynomial coefficients. 
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Fig. 34. — Left panel: Relative apparent speed of trailing components vs. Lorentz factor of 
the main disturbance. Right panel: Apparent speed of trailing components vs. deprojected 
distance from the core. 
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Fig. 35. — Rate of superluminal ejections in the rest frame of the AGN vs. average jet 
Lorentz factor. 
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Fig. 36. — Polarized intensity images of 3C 66A at four successive epochs when highly 
superluminal component B7 is detected. The polarized intensity peak is 19.3 mJy beam -1 . 
The size of the beam is indicated in Table 1. Contour levels increase by a factor of 2, with 
the lowest contour corresponding to 8% of the peak polarized intensity. The line segments 
show the direction of the electric vectors. The solid lines indicate motion of component B7 
and C2 on the total intensity images. 
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Fig. 37. — Polarized flux density of the knots I (solid line) and t (dashed line) in the radio 
galaxy 3C 120 at different distances from the core. The numbers indicate values of the 
EVPA. 
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Fig. 38. — Relative right ascension and declination vs. epoch for components Al (circles) 
and A2 (triangles) in OJ 287. The lines show a straight-line approximation to the temporal 
behavior. 
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Fig. 39. — Total flux density of knots Bl (solid lines) and B2 (dashed lines) in the quasar 
3C 273 at different distances from the core. 



-96- 



15 



10 



a 



5 - 



i r 



i i r 



i r 



B2 



bs t 



3C 273 



I O 1 



Bl 



bl 



b2 





-140 



J I I I I I I I I I I I I L 



-130 -120 

<0>° 



-110 



Fig. 40. — Apparent speed vs projected position angle of superluminal components in the 
quasar 3C 273. 
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Fig. 41. — Total (contours) and polarized (gray scale) intensity images of 3C 273. The total 
intensity peak is 1710 mJy beam -1 , the polarized intensity peak is 90 mJy beam -1 , the beam 
is 0.39x0.14 mas 2 at P.A.=— 6.5°, and contours increase by a factor of a/2, with the lowest 
contour corresponding to 0.25% of the peak. The line segments indicate the electric vector 
direction. The figure can be found at the web-site www.bu.edu/blazars/multi.html. 
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Fig. 42. — Total intensity images with polarization vectors of the quasar PKS 1510—089 from 
1998.94 to 1999.55 (left panel: S peak =1520 mJy beam" 1 , S£ eak =107 mJy beam" 1 ) and from 
1999.76 to 2000.26 (right panel: S pcak =1700 mJy beam" 1 , S£ eai =147 mJy beam" 1 ). The 
beam is 0.40x0.15 mas 2 at P.A.=-5°, contour levels are 0.3, 0.6,. ..76.8, 90% of the total 
intensity peak. The figure can be found at the web-site www.bu.edu/blazars/multi.html. 
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Fig. 43. — Total flux density of knots C9 (asterisks), CIO (circles), Cll (crosses), and C12 
(triangles) in the quasar 3C 345 at different distances from the core. 
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Fig. 44. — Total (contours) and polarized (gray scale) intensity images of 3C 454.3. The 
total intensity peak is 1860 mJy beam -1 , and the polarized intensity peak is 71 mJy beam" 1 . 
The beam is 0.31x0.15 mas 2 at P.A.=— 6°; contour levels increase by a factor of \pl with the 
lowest contour corresponding to 0.25% of the peak. The line segments indicate the polariza- 
tion direction. The figure can be found at the web-site www.bu.edu/blazars/multi.html. 



